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Abstract 

This paper is focused on the generalized Forchheimer flows of slightly compressible fluids in 
porous media. They are reformulated as a degenerate parabolic equation for the pressure. The 
initial boundary value problem is studied with time-dependent Dirichlet boundary data. The 
estimates up to the boundary and for all time are derived for the L°°-norm of the pressure, 
its gradient and time derivative. Large-time estimates are established to be independent of the 
initial data. Particularly, thanks to the special structure of the pressure’s nonlinear equation, 
the global gradient estimates are obtained in a relatively simple way, avoiding complicated 
calculations and a prior requirement of Holder estimates. 


1 Introduction 

In studies of fluid dynamics in porous media, Darcy’s law is ubiquitously used. However this linear 
relation between the velocity and pressure gradient is not always valid. The deviation from Darcy’s 
law is well-known when the Reynolds number increases mm- Such a deviation was even noticed 
in early works by Darcy [3] and Dupuit [6] . Nonlinear alternatives were formulated by Forchheimer 
mm, and were studied extensively afterward in physics and engineering, see [2lll6l[l7l|2ll[26] and 
references therein. In contrast to the vast mathematical research of Darcy’s flows, see e.g. [25], 
existing mathematical papers on Forchheimer flows are much fewer and came much later. Even 
rarer are the ones for compressible fluids. (See, e.g., [Tll^llOj for more introduction to Forchheimer 
flows.) Generalized Forchheimer equations were proposed [HIMIllIs] in order to cover a general 
class of fluid flows in porous media formulated from experiments. In our previous paper |I2] . we 
derive interior estimates for generalized Forchheimer flows of slightly compressible fluids. In this 
article, we focus on the spatially global estimates, i.e., on the entire domain. Since the equation 
for pressure is of degenerate parabolic type, finding L°°-bounds for its gradient, in general, is a 
difficult task and requires much work, see, for e.g., [5l[l8l[l9l[2l| . We will show that the gradient 
estimates for the particular flows under the current study can be obtained in a relatively simple 
way thanks to their equations’ specific structure. 
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In this paper, we consider the initial boundary value problem for the pressure in a bounded 
domain with time-dependent Dirichlet boundary data. Our goal is to obtain bounds for the pres¬ 
sure’s gradient and time derivative in terms of initial and boundary data. Aiming at studying the 
long-term dynamics of the system, we also emphasize the bounds for large time. The standard 
technique from see also mm, requires a Holder estimate first, which, by itself, is not a small 
task. Moreover, the estimates obtained in m for degenerate equations depend on the initial value 
of the gradient. This will result in long-time estimates for the gradient that are dependent on the 
initial data. That means that the long-time dynamics cannot be reduced to a much smaller set, say, 
the global attractor (c.f. [2ni[23]). In this paper, we will demonstrate that these two setbacks can 
be overcome without overcomplicated calculations. First, we extend the techniques for parabolic 
equations in [Hj to our degenerate one in a suitable way. Such extension is possible due to the 
equation’s special structure. Second, by doing analysis on the entire domain, we avoid the spatially 
local consideration, and hence, the requirement for the Holder estimates. Finally, localizing the es¬ 
timates in time and utilizing uniform Gronwall-type inequalities remove the estimates’ dependence 
on the initial data. 

The paper is organized as follows. In section [2l we recall basic facts about generalized Forch- 
heimer equations, and prove a global embedding of Ladyzhenskaya-Uraltseva type. It is different 
from the versions in [12tll4| which are localizations at an interior or a boundary point. This will 
contribute to simpler proofs for the global estimates in section [6l In section [3l we review some 
relevant estimates from our previous works [111112] which will be used repeatedly. In section H] 
we use De Giorgi’s iteration to obtain (spatially) global estimates. Time-local inequality ()4.2p is 
“quasi-homogeneous” in terms of L" ^-norm of the solution. This improves estimates in our previ¬ 
ous works, and is similar to a recent improvement on interior estimates in [22]. Explicit estimates 
in terms of initial and boundary data are established in Theorem 14.41 Section [5] is devoted to 
estimating the maximum of the gradient’s modulus on the boundary. Unlike the results in [14] . 
we are able to localize the estimate in time in Theorem 15.21 hence, it is independent of the initial 
gradient. This is essential to finding the gradient bounds for large time, by combining it with 
certain uniform Gronwall-type estimates. This is demonstrated in Gorollarv 15.31 In section (U we 
establish global L'^-estimates for the gradient for all s > 0 and all positive time. This requires only 
moderate regularity for the initial data, say, L" and with an appropriate a > 0. Section [7] 

contains L°°-estimates for the gradient. Thanks to the special form of the resulting equations (|7.2p 
for the gradient, it is possible to apply the De Giorgi technique, see Theorem 17.11 Section [8] has 
the same results as section [7] but for the time derivative. The main estimates are in Theorems 18.21 
and 18.31 


2 Preliminaries 

Consider a fluid in a porous medium in space M"'. For physics problem n = 3, but here we consider 
any n > 2. Let x E M” and t E M be the spatial and time variables. The fluid flow has velocity 
v{x,t) £ IR”", pressure p(x, t) E M and density p{x,t) £ [0, oo). 

The generalized Forchheimer equations studied in [HleHillllS] are of the the form: 

g,(|t,|)^; = -Vp, (2.1) 

where g{s) > 0 is a function defined on [0, oo). When g{s) = a, a + fls, a + fis + a + 
where a, 171,^^11 are empirical constants, we have Darcy’s law, Forchheimer’s two-term, three- 
term and power laws, respectively. 
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In this paper, we study the case when the function g in ()2.ip is a generalized polynomial with 
positive coefficients, that is. 


g{s) = aos"“ + + ... + oats"^ for s > 0, 


( 2 . 2 ) 


where > 1, the powers ao = 0 < ai < ... < ajsf are fixed real numbers (not necessarily integers), 
the coefficients oq, ai, ..., oat are positive. This function g{s) is referred to as the Forchheimer 
polynomial in equation (|2.ip . 

From (12.1|] one can solve v implicitly in terms of Vp and derives a nonlinear version of Darcy’s 
equation: 

v = -K{\Vp\)Vp, (2.3) 

where the function K : [0, oo) —)> [0, oo) is dehned by 




1 


with s = s(^) > 0 satisfying sg{s) 


for ^ > 0. 


(2.4) 


In addition to m, we have the continuity equation 


i^ + V-{pv) = 0, 


(2.5) 


where number (/> G (0,1) is the constant porosity. Also, for slightly compressible fluids, the equation 
of state is 

^ , with K = const. > 0. (2-6) 

dp K 

From (j2.3p . (j2.5p and (12.6p one derives a scalar equation for the pressure: 


= kV • (iF(|Vp|)Vp) + iF(|Vp|)|Vp|2. (2.7) 

On the right-hand side of ()2.7p , the constant k is very large for most slightly compressible fluids 
in porous media m. hence we neglect its second term and by scaling the time variable, we study 
the reduced equation 

^ = V-(A(|Vp|)Vp). (2.8) 

Note that this reduction is commonly used in engineering. 

Our aim is to study the initial boundary value problem for equation (|2.8I) in a bounded domain. 
Here afterward 17 is a bounded, open, connected subset of M"’, n = 2,3,... with boundary 
F = dU. Throughout, the Forchheimer polynomial g{s) in (|2.2I) is fixed. The following number is 
frequently used in our calculations: 


OiN 

1 -(- CX]\f 


G(0,1). 


(2.9) 


The function K{^) in (12.41) has the following properties (c.f. [Il[9]): it is decreasing in ^ mapping 
^ G [0, oo) onto (0, l/oo], and 


(i + 0“ (1 + 0“’ 

(2.10) 

- 1) < K{^)e < 

(2.11) 

- aA(0 < iF'(0e < 0, 

(2.12) 
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where di, d2, d^ are positive constants depending on g. As in previous works, 
H{^) defined by 

re 

H{^) = / K{y/s)ds for ^ > 0. 

Jo 

It satisfies < H{^) < thus, by (|2.1ip . 

doif-" - 1) < H{0 < 2d2^2-“. 


we use the function 


(2.13) 


The following parabolic Poincare-Sobolev inequalities are needed for our study. For each T > 0, 
denote Qt = U x (0,T). We define a threshold exponent 


an 


a* = 


2 — 0 

Lemma 2.1 (cf. [12], Lemma 2.1). Assume 

a >2 and a > a* 
Let 


p = a 1 + 


2 — 0 


n 


— a. 


Then 


(2.14) 

(2.15) 

(2.16) 
(2.17) 


\\uh.iQ,)<C{l + 5 T)eP[[u]], 

where 6 = 1 in general, 5 = 0 in case u vanishes on the boundary dU, and 

[[n]] = esssup ||n(-,t)||ia(( 7 ) + (" f f |u(x, t)|"“^| Vu(x, . (2-18) 

[o,T] ^JoJu ’ 

The next is a particular embedding with spatial weights from Lemma 2.4 of m (see inequality 
(2.28) with m = 2 there). 


Lemma 2.2 (cf. [13] . Lemma 2.4). Given W{x,t) > 0 on Qt- Let r be a number that satisfies 

2 n 


Set 

Then 


n + 2 


Q = Qir) = 4(1 — 1/r*). 


<r <2. 


(2.19) 

( 2 . 20 ) 


e 2 — r 

+ ess sup ( / W{x,t)~~ Xsupp uix, t)dx) (2.21) 

tefo.Tl ^ Ju ' ’ 


te[o,r] 

where 5 = 1 in general, 6 = 0 in case u vanishes on the boundary dU, and 


[[u]]2,W\T = ess 
[0 


snp \\u{-,t)\\L^u) + ilL Wix,t)\Vu{x,t)\'^dxdty. (2.22) 


Above, supp/ denotes the support set of a function /, and XA denotes the characteristic func¬ 
tions of a set A. 

The following embedding is a global version of Lemma 3.7 from m and Lemma 5.4 on page 
93 in [H]. 
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Lemma 2.3. Suppose w is a function on U that satisfies \\iw\ < M on F. 

Let V = maxIlVwp — M^,0}. For each s > 1, there exists a constant C > 0 depending on s 
such that for any k 

[ K{\Vw\)v^+^dx <C max\w-k\‘^ [ Ki\Vw\)\V^w\‘^v^-^dx + CM^ [ Ki\Vw\)v^-^dx. (2.23) 
Ju Ju Ju 

Proof. Let / = K{\Vw\)v^^^dx. Note that u = 0 on r. First, we see that 

I< [ K{\Vw\)v^\Vw\‘^dx = ^ [ K{\Vw\)v^diwdi{w - k)dx. 

Ju Ju 

By integration by parts, 

n « 

/ = — / di{K{\Vw\)v^diw) ■ (w — k)dx 

1^1 Ju 

= — ^ f (K'{\Vw\)^^^^^^^f^^v‘^diW-{w — k)dx— [ KdVwDv^^Aw ■ {w — k)dx 




\Vw\ 


'U 


— s 


IL « 

/ K{\Vw\)(^v^~^didjwdjw'^ ■ diW ■ {w — k)x{v>o}dx- 

i,j=i Ju 

From this and (|2.12l) . it follows that 

I<a [ K{\Vw\)\V‘^w\v^\w — k\dx + [ K{\Vw\)v^\Aw\\w — k\dx 


lu 


lu 


+ s I ^\'V‘^w\\'Vw\'"‘\w — k\x{v>o}dx 


lu 


<C K{\Vw\)v^\V^w\\w — k\dx + C / iF(|Vt(;|)|V"^ry|?;* ''{v + M‘^)\w — k\dx. 


lu 


lu 


Hence, 


I<C I K{\Vw\)v^\V‘^w\\w-k\dx + CM'^ I K{\Vw\)\V^w\v^-^\w - k\dx 




lu 


lu 


= C K{\Vw\)v 2 • |V^rc|t; \w— 


lu 


\w — k\dx + C f KdVwDv^^ ^\V‘^w\\w — k\ ■ M‘^dx. 

Ju 


This last inequality and Cauchy’s inequality imply that 

I<J K{\Vw\){^^ + C\V‘^w\‘^v^-^\w-k\‘^)dx + C [ K{\Vw\)v^-\\V^w\‘^\w-k\‘^+ M^)dx 


lu 


= -I 
2 


+ C [ K{\Vw\)v^-^\V^w\‘^\w - k\‘^dx + CM^ [ K{\Vw\)v^-^dx. 
Ju Ju 

Therefore, we obtain (I2.23F 


□ 


The following is a generalization of the convergence of fast decay geometry sequences in Lemma 
5.6, Chapter II of [14]. It will be used in the De Giorgi iterations. 
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Lemma 2.4 (cf. [13] . Lemma A.2). Let be a sequence of non-negative numbers satisfying 

m 

k=l 

where B > 1, > t) and > 0 for k = 1,2,... ,m. Let g, = mln{gLk '■ 1 < k < m}. If 

Yq < ■ 1 < k < m} then llmi^ooYi = 0. 


3 Previous results 


We study the following initial boundary value problem (IBVP) for p{x,t): 

'^ = V-(A:(|Vp|)Vp) inC/x(0,oo), 

^ p{x,0) = pq{x) in U, 

p{x,t) = ipixjt) on r X (0, oo). 


In order to deal with the non-homogeneous boundary condition, the data ipix, t) with x E L and 
t > 0 is extended to a function 'I>{x,t) with x G U and t > 0. Throughout, our results are stated 
in terms of 'h instead of iJj. Nonetheless, corresponding results in terms of fj can be retrieved as 
performed in [9]. The function T is always assumed to have adequate regularities for all calculations 
in this paper. 

Solutions. It is proved in section 3 of m that m possesses a weak solution p{x, t) for 
all t > 0. It, in fact, has more regularity in spatial and time variables, see [4]. For the current 
study, we assume that solution p{x, t) has sufficient regularities both in x and t variables such 
that our calculations hereafter can be performed legitimately. Specifically, we assume p,^p,pt E 
C{U X (0, oo)), the Hessian matrix of second spatial derivatives E C{U x (0, oo)), and the 
function t —)• p{-,t) is continuous from [0, oo) to n L°‘(U), for an appropriate a > 0 

which will be determined for each type of estimate. Some particular estimates may require much less 
regularity such as the L°°-estimates in section [H The obtained estimates also hold without the C^- 
requirement by using an approximation process, see m- However, to avoid further complications, 
we do not perform such approximation here. 

Generic constants notation. Hereafter, the symbols C and C are used to denote positive 
numbers independent of the initial and boundary data, and the time variables t, Tq, T; it may 
depend on many parameters, namely, exponents and coefficients of polynomial g, the spatial di¬ 
mension n and domain U, other involved exponents a, s, etc. in calculations. The values of C and 
C may vary from place to place, even on the same line. 

Functions and Norms. Throughout, whenever unspecified, the norm Lebesgue or Sobolev 
norms mean for the whole domain U. Also, for a function /(x, t), we use t{t) to denote the functions 
t f{-,t)- For example, ||p(t)||L 2 = \\p{-,t)\\LHu): max^^r] IIp(0IIl°° = max^gp^r] \\p{-,t)\\L^(u)- 

We recall some relevant results from m- For a > 1, we define 


A{a, t) 



|V4/(x,t)| 


a{2 —a) 

2 dx 


2(cy — a) 
q( 2 —a) 





4't(x,t)|"dx 


a —a 
q( 1—a) 


(3.2) 


A(a) = limsupA(a, t) 

t—>oo 


I3{a) = limsup[A'(a, t)] . 

t—^OO 


for t > 0, and 


and 


(3.3) 
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Also, define for a > 0 the number 

a = max {a, 2, a* }. (3-4) 

Whenever /3{a) is in use, it is understood that the function t A{a,t) belongs to C'^((0,oo)). 
For a function / : [0, oo) ^ M, we denote by Envf a continuous and increasing function 
F : [0, oo) —)• M such that F{t) > f{t) for all f > 0. 

Let p{x,t) be a solution to IBVP (13.Ih . and denote p = p — 

Theorem 3.1 (cf. [11], Theorem 4.3). Let a > 0. 

(i) For all t > 0, 


'U 


x,t)|"(ix < Cl 1 + / \p{x,0)\'^dx + [EnvA{d,t)]°‘-<^ ). 


lU 


(ii) If A(d) < oo then 


limsup / |p(x,t)|“dx < C(l + A(a)s-a). 
*—>■00 Ju 


(iii) If /3(a) < oo then there is T > 0 such that 

/ \p{x,t)\°‘dx < C(l + /3(a)s-2“ + A{d,t)^^) forallt>T. 

Ju 

For gradient and time derivative estimates, we denote 

Gi(t)= [ |V^(x,t)|2dx+ [ [ \^t{x,t)rdx 


(3.5) 


(3.6) 


(3.7) 


lu 


'u 


ro(l-a) 


+ 


|4't(x, t)Y°dx 


'-JU 


ro 


G2{t) = / \VMx,t)\^dx+ / \Mx,t)\‘^dx, G3it) = Giit) + G2{t). 


'u 


lu 


with ro = ^2-2){n+i)-n • i > 0, recall from (4.20) in [TT] and from (3.25) in [9] that 

f f H{\Vp\)dxdT < C f p^(x,0)dx + C f Gi(r)fir, (3.8) 

Jo Ju Ju Jo 

[ II{\Vp\){x,t)dx+ [ I \pt{x,T)\'^dxdT < I [H{\Vp{x,0)\) + p‘^{x,0)]dx + C [ G^{T)d‘ 

Ju Jo Ju Ju Jo 


T. 

(3.9) 


Below are estimates for gradient and time derivative when time is large. 

Theorem 3.2 (cf. [T2|, Corollary 3.3). Let a > 2. 

(i) For t > 1, 

f f H{\Vp{x,T)\)dxdT < C(l f \pQ{x)\°‘dx [EnvA{a,t)]°‘-°- + f Gi(r)drV (3.10) 
Jt-l Ju ^ Ju Jt-1 7 

[ [ Pt{x,T)dxdT < C(l-I f \po{x)\°‘dx -I [EnvA{a,t)]°‘~°' + [ Gs{T)dT\ (3.11) 

Jt-l/2 Ju ^ Ju Jt-l ' 

(ii) If A{a) < oo then 

limsup f f 77(|Vp(x,r)|)(3xdr < Cf^l + A(a)“^ + limsup f Gi{T)dT\ (3-12) 

*^oo Jt—l Ju ^ *—>00 Jt—1 ' 
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limsup / / pf{x,T)dxdT<c(l + A{a)°‘-<^+limsiip / G^{T)dT\. 

t^oo Jt-1/2 Ju ^ t^QO Jt-1 ' 

(iii) // /3(a) < oo then there is T > 0 such that one has for all t>T that 

f f H(\Vp(x,T)\)dxdT <c(l + (3(a)°‘-'^°- + A{a,t — 1 )^^ + f Gi(T)dT\ 
Jt-i Ju ^ Jt-i ^ 


(3.13) 


f 

lt- 1/2 


J p^(x,T)dxdT < c(^l + f3(a)°‘-'^°'+ A(a,t — + J G^(T)dT 


(3.14) 

(3.15) 


4 -estimates 


In this section we estimate the L°°-norm of the pressure. We will focus on estimates for p instead 
of p. First, we give a local (in time) L°°-estimate which does not depend on the initial data’s 
L°°-norm. 

Let p{x, t) be a solution to IBVP (|3.ip with given data po{x) and f^ix, t). Let p = p — ^, then 
it satisfies 

dp 


= V • iK(\Vp\)Vp) inUx (0, oo), 

\^p{x,t) = 0 on r X (0, oo). 

Theorem 4.1. Let a satisfy (12.151) . /f Tq > 0, T > 0 and 9 E (0,1) then 


(4.1) 


sup 

[To+0T,To+r] 




{m iipiil<.(,ix(to,to+t)) 


L°‘{Ux{To,To+T)) 


+ 


(9T) ||p||/a(f/x(To,ro+T)) + r2(“ “)(i+' 52) ||p||(^ x(Tq.Tq+T)) 


+ {T^\\V^\\L‘^<i^uxiTo,To+T)) +^~ll^t|lL“U(!7x(ro,To+T)))^+^3 l|p|lL“(r/x(To,To+T)) 


«4 


+ (r2a ||V4'||ia,i(^x(ro,To+T)) “ ll^t||L“U(r/x(To,To+T)))^^^^^^^^^l|p|lLt(l/x(To,To+T))}- 

(4.2) 


where ri = q;( 1 + (2 — o)/n) — a, numbers pi and qi are conjugates of each other that satisfy 

1 < Pi < or equivalently, 1 H—-—;- — < qi = ——— < oo, 

a(l/a* — 1/a) Pi — 1 

and hi, 82 , hs, 84 are positive numbers defined by 

8 i = l-a/ri, 82 = al(a - a) - afri, 

<53 = 1/Pi - a/ri, 84 = a/((a - a)pi) -ajrx, 

Proof. Without loss of generality, we assume Tq = 0. Let A: > 0, define = max{p — /CjO}, and 
denote by Xk the characteristic function on the set supp p^^'i - the support set of p^^\ 

Let C = C(t) be a smooth function on M satisfying 0 < C < 1 and Ct ^ 0 oa M, and </ = 0 on 
(- 00 , 0 ]. 
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Multiplying the partial differential equation (PDE) in (13.111 by integrating over U 

and using integration by parts, we have 


1 

« Ju 


[ + (a - 1) [ K{\Vp\){Vp-Vp^^'^)\p^^^\‘^-^Cdx = - [ ^t\P^^'^\'^-^Cdx. 

Ju Ot Jjj Jjj 


For the second integral on the left-hand side, 

X(|Vp|)(Vp • = K{\Vp\){Vp ■ + i^(|Vp|)(V^ • 

Note that Vp-Vp^^^ = and jV'!' • Then 

K(|Vp|)(Vp- > K(|Vp|)|#)|“-2(|VpW|2 _ i|VT|2 - 

= lKi\Vp\)\Vp^^^\V^r-^ - iK(|Vp|)|VT|2|;pW|“-2. 


Therefore, 


a dt 

; 


Q|p(fc)|a a—1 

(dx H--— 


lu 


K(|Vp|)|Vp(^)|2|pW|“-2cdx 


'u 


K(|Vp|)|VT|2|pW|“-2(dx+ / 


(4.4) 


/[/ 


<C [ |VT|2|;p(^)|“-2(dx+ [ \^t\\P^^^r~\dx. 

Ju Ju 


The last inequality uses the fact that K{^) is bounded. 

For the second integral on the left-hand side again, we have from (|2.10ll and triangle inequality 
that 

12 di|VpW|2 di|Vp(^)|2 


K{\Vp\)\VpJ''^\^ > > 

Uvpuivp I - 

(di/2)(|Vp(*=)| + |VT|)2 - fii|V^|2 


(l + |VpW| + |VT|)“ 


> 


> 


(l + |VpW| + |VT|) 
di (|Vp(^)| + |VT|)2 


2 (1 + |VpW| + |VT|)“ 
It is elementary to show for ^ > 0 that 

. 1 


-di|VT| 


2—a 


> ^(,c2-a _ 


(l + 0“ “ 2“ 


Thus, 


A-(|Vp|)|Vp'‘>P > ^({|ViiW| + |VW|f-“ - 1) - d,|V4'p-“ 

Therefore, it follows this, (j4.4h . and Young’s inequality that 

/ / |Vp('=)|2-“|p(^)|"-2Cdx 

Ju Ot Jjj 

<C [ (1 +|V^|2-“ +|VT|2)|pW|“-2(dx + (7 [ |Tt||p(^)|“-iCdx 

Ju Ju 

<C [ (1 +|V^|2)|pW|“-2(dx + (7 [ \^t\\P^'^'>\^-\dx. 

Ju Ju 
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Applying the product rule to the first term, we have 

4 - [ \P^^^\'"Cdx+ [ 


dt 


lu 


lu 


<C \p^’^^\^Ctdx + C / ((l + |V^|2)|pW|“-2c + |^i||pW|“-V)dx. 


'u 


lu 


Applying Young’s inequality to the integrand of the last integral yields 

4- f < C [ l^^rCtdx 

dt 


'u 


lu 


lu 


+ £ 


[ i#)rcdx+c£i-“/2 f (I+ \s/^>\»)xkCdxdt+c£^-^ [ i^trxkCdxdt. 

Ju Ju Ju 


Integrating form 0 to t for t G [0, T] and taking the supremum give 


fl 

lo Ju 


sup [ \p^^'>\‘^Cdx+ [ [ < C 

[o,r] Ju Jo Ju 

+ eTsup [ |p(^)rC(^x + Ce ^-“/2 f f (i + |v^|“)XfcC(^xdt + Ce^-“ [ [ l^t^XkCdxdt. 

fO.Tl Ju Jo Ju Jo Ju 


[0,T]JU 

Selecting e = 1/(2T), we obtain 

rT 


\p^^'^\‘^Ctdxdt 
T 

0 JU 


r f 

0 Ju 


sup [ |p(*’)|“C(^x+ [ [ < C 

[ 0 ,T] Ju Jo Ju 

+ Cr“/2-i r [ XkCdxdt + r I \V^\^XkCdxdt 

Jo Ju Jo Ju 

Applying Holder’s inequality to the last two integrals gives 

sup [ |pW|“C(^x+ [ [ |VpW|2-“|p(*’)|“-2Cdxdt < C 

[o,r] Ju Jo Ju 

+ r [ XkCdxdt 

Jo Ju 


\p^’^^\‘^Ctdxdt 

rT 


+ CT 


iQ:—1 


'0 Ju 


\^trXkC,dxdt. 


0 Ju 


\p^^'^\^Ctdxdt 


(4.5) 


+ C{T-/^-^\\\V^rhn(u><(0.T))+T^-^mr' 


|L'?1(C/X(0,T)) -I-J lll'^tl \\Lii{Ux(0,T))) Hi 

Let Mq > 0 be fixed which will be determined later. For f > 0, define 

ki = Mo{l-2-^), L = 0T(l-2-*). 


r s 1/pj 

XkCdxdt ] 

0 JU 


(4.6) 


Then to = 0 < ti < ... < 6 T and limj_^oo L = 6 T ; feo = 0 < /ci < A :2 < • • • and limj_>.oo ki = Mq. 

For z,j > 0, we denote 

Qi = Ux{ti,T) and Aij = {{x,t) ^ Qj ■ p{x,t) > h}, Ai = Ai^i. (4.7) 

For each f, we use a cut-off function Ci{t) with = 1 in [tj+i,T] and Ci = 0 on and 

C C'2*+i 


l(0)tl< 


L-l-l ti 


er 


(4.8) 
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for some C > 0. Applying (14.5p with k = ki+i and C = Ci gives 

sup [ \p^^^+^^\°Cidx+ [ [ 

fO.Tl Ju Jo Ju 


< 


Jo Ju Jo Ju 


idxdt 


(4.9) 


'0 JU 

+ c(r"/2-i||v4' 


IL“m([/x(0,T)) 


+ r“ ^||^t|lL“M(c/x(o,r))) 


^ \ 1/pl 

/ Xki+iCidxdt] 

0 Ju 


Define 


Let 


Fi= sup [ 

[ti+uT]Ju Jti Ju 


fl=r“/2-l and ^2 = r“/2-l||V4/||2..i(t;,(0,T)) 
Then (14.9p yields 


+ r"“^||'L "" 


Fi < 


rU+i 

Iti JU 


ti+i Ju 


Xki+idxdt + CS2 


Xki+idxdt 


ti+i Ju 


+ CTi 

< +CTi|A,+i,| +CT2 |A,+i,|1/^^ 

Since \\p^’^J\\L^{Ai) > llL-(yii+i,i) > ih+i - fc*)!A^+pilthus 

\A,+iA < {h+i - 
This and (j4.10p imply 


t|lL“9i(!7x(0,T))- 

1/pi 


i 


(4.10) 


(4.11) 


< C2' 


“{((0r)-i+TiMo 




Note that ri is the exponent defined in (12.161) . By Lemma l2.11 

< c(Fy“+Fy<‘>-“>). 

Holder’s inequality gives 

Combining (14.140 with (14.lip . (|4.12p and (I4.13P yields 

1/q 




‘{((OF)-' +£iA4-“) ||pl‘->||t.(A,) +£l''“A/„“"''’'‘||p'‘-'lli2’,‘ 



II “/Pl 

llL“(£li) 

II “/pi 

llL“(Ai) 

}. (4.12) 


(4.13) 


1 l/a—1/ri 

M 


(4.14) 

^||p(fci)|| 

i/pi 


+ 


-1 


l/(Q-a) 


+ j ' ' WP^^J + ^V(«-a)^-«/((a-a)pi) II -(fc,) ||a/j(a-a)pD I 

-l+a/ri|| (fc,)||l-a/ri 

IIP IIl2(^,) , 

where B = 2^. Let Yi = \\p^^J\\L<^(Ai)- Then 

Ti+i < CB^ + ^4^/+'^") , 
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where 




ZI 2 =^2 M, 


1/" A^-l/Pl . 


0 




0 


= £ 2 ' 




--5i 


(4.15) 

(4.16) 


ZI 4 = ^i/(«-)^-«/((«-'^)pi) . 

Take Mq sufficiently large such that 

or equivalently, 

Dj<CY-'\ j = l,2,3,4. 

Then by Lemma I^ITI limj_).oo li = 0, consequently, \p^^°^\'^dxdt = 0, that is 

p{x,t) < Mq in [7 X [6T,T]. 

Repeating the argument above for —p,—V’ instead of p, we obtain 

\P{xH)\<Mo inUx[eT,T]. 

It remains to determine Mq. Since Yq < \\p\\L‘^(ux{o,T))y we have sufficient conditions for (j4.16l) : 

(9r)-=i;A4-* < c||p||2.Tx(0,t)). < c||p||zJ,„,|„_^„, 


(4,17) 


1/" ^ r'll4^l|-'52 (.a-a)p-^ ^ r* 11 I “<54 


£ 2 ' M, 


^ £'\\p\\l°‘\ux{0,T))^ ^2 

Solving these inequalities gives 


Mq > C max 


{m 




^ £'\\p\\lHUx(Q,T))- 


L“(f/x( 0 ,r))) \\p\\l°‘(Ux{0,T))’ ll7'llL“(i7x(0,r))’ 

S^(at-a) 1 S 2 


7"2(c«+5i(c,-a)) IIztII “+'5i(“-“) 

^ ^ \\P\\l‘^(Ux{0,T))^ ^2 


.Q(l/pi+il) II -|| l/pi+^l 


lL“(!7x(0,T))’ ^2 ll■f^llL“({/x(0,^))/ 


_ n _ i54(Q-a)pi 

Q; + (5i(ct-a)p4 II -|| a + ii(Q;-a)p4 ( 


(4.18) 


with an appropriate positive constant C. Choosing Mq to be the right-hand side of (j4.18p with the 
sum replacing the maximum, we obtain from (|4.17p that 


“ ^ I-II i+ii 


sup Mt)\\L^(u) <cUeT) ||p||L-(t/x( 0 ,r)) + r 2 (Wi)“ I 12.11 lo.>ux(q 0 +t)) 
{eT,T] 4 V . 


+ {eT) 


£3 


c-2 . 

I 'r2(c + SPa-a)) ||^|| 0‘ + SP°‘-o^) 

^ 1 '• \\P\\l°‘(Ux{0,T)) 


lL“([/x(0,T)) 

Rewriting the powers in (I4.19p . noticing that 


L“({/x(0,T)) 

,, 1 _ £2 _ „ j 4 (a-a)pi 

(d'”)"7swiiPii;£y?,„.r«+(T”)"'-+‘“"->iiPiiss;:7„t!>|, (“-ig) 


1 I- . « 

-h5i = 1-7 53, - 

Pi a — a 


-7 5i = 1 -7 ^2, 


a 


{a - a)pi 


lL“([/x(o,r))/- 

“7 5i = 1 -7 ^4, 


and also 


\ / (y^ Ck — 2 Q: — 1 

£■ 2 ! < l|V'I'llL“U(r/x(To,To+T)) + ^ “ ll'I't||L“‘Ji(;7x(To,To-l-T))), 

we obtain (14.21) . The proof is complete. 


□ 
































Global Estimates for Generalized Forchheimer Flows 


13 


Remark 4.2. In (j4.2p . the norm \\p\\(ux[To+eT,To+T]) estimated by a sum of homogeneous 
terms in ||j>||L“(r/x(To,To+T)) • Therefore it is appropriate for both small and large ||p||Lc.(( 7 x(ro,ro+r)) • 
This improves our previous interior versions in Theorem 4-1 of \12^ and Proposition 3.2 of ]13^ . 
where the bounds contain an additive constant, which make them more suitable for large value 
\\p\\l°‘(Ux{To,To+t)- Tor doubly nonlinear parabolic equations of p-Laplacian type, the interior esti¬ 
mate of this kind was obtained in l2^ using Moser’s iteration. 

To simplify our future estimates, we use the following weaker version of Theorem 14.11 
Corollary 4.3. Let a, Tq, T, 6 , ri, pi, qi, 6 i, 62 , 63 , ^4 be as in Theorem \4-l\ Then 


sup < c(i + (0r)-vh 

[ro+0T,ro+r] ^ ^ 


1 + ||V'I'||^a9i((7x(ro,To+T)) + \\’^t\\L°^n{Ux{To,To+T)) ^ 

\ 2:3 

1 + l|p||L“(r/x(To,To+r)) j , (4.20) 


where 


a — 1 


a 


zi = 


{a - a)(l + < 54 ) 


, ^2 = 


{a - a)(l + 54 )' 


2:3 = max 


I ’ 1 + <53 / 


(4.21) 


Proof. We apply Young’s inequality in (j4.2jl to the terms involving (OT) T, ||VT|| 2 ,c« 9 i, , 

\\P\\l°‘- 

First, we note that ^2 > <5i > (^3 and 62 > 64 > 63 . 

• For the power of {9T)~^, we have l/Ji>l/((a — a)Ji). 

• The largest power of T is 


a — 2 


a — 2 


a — 1 


a — 1 


Note that 


zi = max • 


a — 2 


2 a(l + 5i) ’ 2 (a — a)(l + ^ 2 ) ’ a(l + ^ 3 ) ’ (a - a)(l + ^ 4 ). 


< 


a — 2 


0 — 2 


< 


0—1 


2 o(l + Ji) 2(0 + (5i(o — a)) 2 (a — a)(l-h 62 ) (o — o)(l + ^ 4 )' 


0 — 1 


0 — 1 


< 


0 — 1 


0 — 1 


0(1 + 53 ) a(l/pi+5i) (a - a)( (^_°)p^ +5i) (o-a)(l + 54) 


Hence zi is as in (|4.21l) . 

• The largest power of || and is 


Z 2 = max 


(r 


1 


a 


+ ^3 (o — o)(l + 54 ) 


}■ 


Since ^^ < 7 -w —% -= ?- T, , ^ . we have zn as in (|4.2ip. 

1 ++ l/pi+<5l (a-a)( (^-^)p^ ++) (a-a)(l+<54) ’ ^ 'I - 

• The largest power of IIpIIl^ is 


_ [ 5i 02 03 03 04 "I 

’TT^’TT^%’TT^’TT^J' 


Note that < 1- Hence 23 is as in (I4.2ip . Therefore, we obtain (14.201) from 

(i+^ . 112 4 ^ 
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We then derive the L°°-estimates in terms of the problem’s data. 


Theorem 4.4. Assume a satisfies (I2.15|] . Let qi, (5i be defined as in Theorem ]^. 1\ and zi, Z 2 , z^ 
be defined as in Corollary \4.3\ 

(i) If 0 < t < 3 then 

\\P{'t)\\L°°(u) < Ct {l+\\po\\lcy(^^) + [EnvA{a, t)] )^’^(l+||V'I'||2,a9i((7x(o,t)) + ll'^’t||L“9i(!7x(o,t)))^^- 

(4.22) 

If t > 1 then 

\\Pi't)\\L°°{u) < C{l + \\po\\Lai^u) + [EnvA{a, t)] )^’^(l + ||V'I'||2,a9i([/x(t-i,t)) + ll'I't||L“'Ji(t/x(t-i,t)))'^^- 

(4.23) 

(ii) If A{a) < oo then 


limsup||p(t)||Loo < C'(l+A(a)“-“)^3(l+limsup(||V4'(t)||i«,i([;x(t-i,t)) + ll^t(i)lli,“9i(t/x(t-i,t))))^"- 

t^coijj) t^OO 

(4.24) 

(iii) If /3(a) < oo then there is T > 0 such that 


\l°°(u) < C(l+(3{a)°‘ 2“ + ||^(a, ^^)^3(l+||V'I'||£,o,5i([7x(t-i,t)) + ll'I't||L“9i(c/x(t-i,t)))^^ 

(4.25) 

for all t >T. 


Proof, (i) Note that a = a. Let t E (0, 3]. Applying (|4.20p for T = t and 0 = 1/2, we have 

\\p{t)\\L°° < C{1 + t )(1 + ||V4'||iagi(^x(o,L) + W'^tWi’^n/ux/op/Y^i^ + llp|lL“(r/x(o,L))^^- (4-26) 

Using (13.5h to estimate ||p||L“([/x(o,t)) iii (14.26p . we obtain (I4.22h . 

For t > 1, applying p4.20p with Tq = t — 1, T = 1 and 9 = 1/2 we obtain 


||p(i)||L°° < C(1 + ||V'I'||iQi9i((7x(t-l,t)) + ||'I't||i,°i9i((7xp-l,t)))^^(l + l|p||L“(C/x(t-l,t)))^^- (4.27) 

Again using p3.5p with noticing that the function EnvA(a,t) increasing in t, then we obtain p4.23p 
from (|4.27p . 

(ii) From (I4.27P we have 


limsup ||p(t)||Loo < ^(l+ limsup(||4't||ic.9i([/x(t-i,t)) + ll^'tII l^w(( 7x(t-i,t))))^^ 
t—^oo t—>-oo 

• (1 + limsup ||p||L«(c/x(t-i,t)))^®- 

t^OO 


By dM]), 


limsup t)) ^ sup / |p(x,t)|"dx < ^(l + A(a)) 

t^oo t^oo Ju 


OL/[OL — a) 


(4.28) 


Thus (I4.24P follows. 

(iii) Using ()3.7I) we have for large t that 


f f \p(x,T)\°‘dxdT < c(l + 13(a) + f A(a,T)°‘-°-dT]. 

It-i Ju ^ Jt-i 


Therefore, p4.25l) follows from p4.27l) and p4.29l) . 


(4.29) 

□ 
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5 Gradient estimates on the boundary 

In this section, we estimate the maximum of |Vp| on the boundary T x (0, oo). It wiii be used 
in estimates of L^-norms for the gradient in section [ 6 ] and its L°°-norm in section [71 We extend 
Ladyzhenskaya-Uraitseva’s technique m in order to have better estimates for iarge time. 

We rewrite the PDE in (13.ip for p in the non-divergence form as 


Pt ^ ^ -^ijPxiXj — 0, 


(5.1) 


where 


Aij = Aij{x,t) = K{\Vp{x,t)\)5ij + ^ ^ ^\ xi{x, t)px^ (x, t), i,j = 1 , 2 , ...,n. 

Thanks to (|2.10p ~ (l2.12p we hnd that 

I Ail < (l + a)i^(|Vp|), 

n 

(1 - a)K{\Vp\)\y\^ < Y, Aijywj < K{\Vp\)\y\^ Vy G M" 

*j=i 

By (12.101) . we have more expiicit reiations: 

lAjI < (1 + a)d2(l + |Vp|) 


(5.2) 

(5.3) 


(1 - a)di(l + |Vp|)-“|yp < Aijyiyj < ^ 2(1 + |Vp|)-“|y| 2 , Vy € M" 

*j=i 

We wiii estabiish boundary estimates in case the boundary is flat first. For generai boundary, 
we wiii flatten it out and hence transform equation (15.11) to a different, but simiiar, one. To prepare 
for this transformation, we consider a more generai PDE: 


def 


f^p=pt- Y AijPxiXj + Y 
*d=i 


2=1 


where Aj = Aij{x, t), bi = bi{x, t) for = 1 , 2 ,..., n. 
Define for ^ > 0 the function 

^(0 = 


1 


(1 + 0 “' 


Simiiar to (15.20 and (15.31) . we assume that Aj satisfy 

n 

(Y |A,P)'/2<ciA(|Vp|), 


*.i=i 


C2A(|Vp|)|yp < Y AijyiVj < C3A(|Vp|)|y|^ Vy G 
*,j=i 

where ci, 02,03 are positive constants. In addition, we aiso assume that 


(5.4) 


(5.5) 

(5.6) 


\b{x, t)| < C 4 A(|Vy|) with b = {bi, ■ ■ ■ ,br, 


(5.7) 
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for some positive constant C 4 . Note, particularly when y = e^, a unit vector in the standard basis 
of M"", we have from (15.61) that 

Afc > C 2 i?(|Vp|), k = l,2,...,n. (5.8) 

Denote the ball in R” centered at the origin with radius R by Bji, and the upper half ball by 
B^, i.e., B^ = Br n {x = {x',Xn) '■ x' G R”“^, Xn > 0}. Also, denote Fr = dB^ n {xn = 0}-the 
flat portion of the boundary of B^. 

We start with estimates on a flat boundary. 

Proposition 5.1. Let T > 0 and R > 0. Let p{x, t) G C^’^{B'I^ x [0, T]) n C‘^’^{B^ x (0, T]) be a 
solution of (j5.4p in B^ X (0, r], with Aij and bi satisfying (l5.5P " (j5.7p in B^ x [0,T]. Then 

max \Vp{x,t)\ < C{1 + R~‘^)e^ ex.p (C' max |p|) 

rfl /2 X [eT,T] B+ X [eT/2,T] 

•(1+ _ max + _max |V4'|^ + _ max |V^'I'|), (5.9) 

B+ X [0,T] B+ X [o,r] B+ X [o,r] 

for any 9 G (0,1), where C, C are positive numbers independent of R, T, 6, and 

ho = (5.10) 

2 — a 

Proof, (a) Let v = —1 + e'^°^ with the constant kq > 0 chosen later. Then 

u = 0 on dB^ n {xn = 0}. (5-11) 


For i, j = 1, 2,... , n, we have 


vt = K.Qe'^^Ppt, Vxi = Koe'^°Ppxi, VxiXj = Koe''^°PpxiXj + ^le'^^^pxiPxy (5.12) 

Substituting (j5.12p into ()5.4p gives 

n 

+ Tt - ^ Aij^KQ^e~'"°PvxiXj + '^xiXj - KopxiPxj^ + b ■ (kq + VT) = 0. 


*J=i 


Multiplying by Koe^°P, it shows that v{x,t) solves the equation 

^2 KQp- - 


Vt - AijVxiXj + 5 • Vu = - ^ AijKle'^°^pxiPxj + Koe'"°^( ^ Aij'itxiXj - T* - 6 • VT). 


i,j=l *,j=l 

Thus, by (|5.6p and Cauchy-Schwarz inequality, we have 


hi=i 


n n 

£v < -C2Ky°PK{\Vp\)\Vp\^ + Koe^°p[{Y, IA,f)^/'( |T,,,,|2)i/2 ^ ^ 

i,j=l *,i=i 

Using inequality |Vpp = |Vp — VT P > |Vpp/2 — |VTp, and inequalities (15.51) . (|5.7I) . we get 
Cv < -yKge''°PJF(|Vp|)|Vpp+ c2K^e^°PiF(|Vp|)|VT|2 
+ Koe'‘°p[.A(|Vp|)(ci|V2T| + C4|VT|) + jTi] 
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Let 


Ml = max (ci|V^'I'| + C4|V'I'|), M 2 = max M 3 = 02 max 

B+x[0,T] B+x[0,T] B+x[0,T] 

Then 

Cv < -^Kle^°PK{\Vp\)\Vp\'^ + ko{M 3 Ko + Mi)e^°PK{\Vp\) + Koe^^PM 2 . (5.13) 
Let e > 0. We write M 2 = K^Vpl) ■ M2(l + |V|?|)“ and estimate 

M2(1 + |Vp|)“ < M2(1 + \Vpn <M 2 + + e\Vp\^ 

From this and (15.13^ . we find that 

Cv < -(^Ko -e)ACoe"°^'i^(|Vp|)|Vp|2 

2 / (5.14) 

+ Koe'^°Pk{\Vp\){KoM 3 + M 1 + M 2 + S-‘^C 2 -a)MV( 2 -a)y 

(b) We localize the above calculations. Let C = C(^) ^) ^ [Oj 1] be a cut-off function on Br x (0, T) 
with C = 0 on {Bji x [0T/2,T]) U {dBji x [0,T]), and Cix,t) = 1 on Bii ^2 ^ {G'C,T), and satisfy 

|vcl<f, |WCI<§. IC,I<§. 

where Cq > 0 is independent of R, 9, and T. Denote by x = x(i) the characteristic function of 
[9T/T, T]. 

Let w = Then 

Cw = CCv + vC{C) - 2 CiAVv) ■ VC. (5.15) 

For the last term, by (15. 5 p and (|5.12l) : 

|2C(^Vn) • VCI < 2ciC^(|Vp|)|Vn||VC| < CiCKoe"°PiL(|Vp|)|Vp| 

< Koe"°P.^(|Vp|) [CiClVpIx + CiClVTl], 

where Ci = 2ci. Thus 

|2C(^Vn)-VC| <ACoe"°^^(|Vp|)(eC"|Vpp + e-ic2x + CiC|VT|). (5.16) 

For the second term on the right-hand side of (|5.15l) . 

n 

\vC{e)\ < |u|2C|G| + \v\iY, IAill(C')x,x,| + 2C|5||VC|) < |n|2C|Ci| +C'2|n|xi^(|Vp|), (5.17) 

*J=i 

where C 2 = 6*3(1 -|- R~'^), with 6*3 > 0 independent of R, 9, T. Treating the term 2|n|C|Ct| in (|5.17p 
the same way as we did for M 2 in (j5.13p : 

2|n|C|Ct| < 64 |n|C < C'4|n|CiL(|Vp|)(l + < k{\Vp\){CMC + C4|u|C|Vpr), 

where 6*4 = 2C'o(0T)“^. Note that 

C4|n|C|Vp|“ = C4|n|C^““Noe^°^]”“/^ • e^°P/^C,\Vp\Y. 
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Applying Young’s inequality with power 2/(2 — a) and 2/a yields 

2|u|C|Ct| < i^(|Vp|)[c4|u|C+(C4|u|C^-“)2/(2-'^)(eAtoe"°n-“/("-“)] +eKoe^°Pk{\Vp\)C^\Vp\^ (5.18) 
It follows from (I5.17P and ()5.18p that 
\v\£{e) < eKoe^°Pk{\Vp\)C^\Vp\^ 

+ k{\\/p\) [C4|u|C + + C2\v\x] ■ (5.19) 

Using (15.141) . (I5.19P and (|5.16l) in (I5.15p . we obtain 
Cw < -C"(|«o - 3 e)Koe^°Pk{\Vp\)\Vp\^ 

+ Koe^°Pk{\Vp\) [c\koM3 + M 1 + M 2 + + CiC|V^|] (5-20) 

+ ^(|Vp|)[C2x|u| + C4|u|x + (C4|u|x)'/('-“)(£ACo)-“/('-“)e"°IPl“/(2-a)]. 

Choose £ = 1 and kq = 8/02- Then 

Cw < Koxe^^^^^k{\Vp\) [kqMs + M1 + M2 + + Cf + Ci|VT| 

+ X(|Vp|)[C2xbl + CMx + 

Let M4 = max^+^T] 1 ^ 1 - that Define 

Ms = Koe''o^4(KoM3 + Ml + C0M2 + + Cf + Ci max |VT|) 

B+x[0,T] 

^ g«oM 4 { 2 +a)/( 2 -a)^^ ^ ^-a/( 2 -a) ^ + 04 + C|^“). 

We obtain Cw < MsiLdVp]). Note that 


(5.21) 


Ms < (ko^s + Ml + C 0 M 2 + + Cf + Ci max |VT|) 

B+x[0,T] 

^ gK0M4(2+a)/(2-a) ^ Kq “/^^■“^)(1 + ^2)2(1 + ) < Mg, 


where Mg is explicitly defined by 


Afg = C'(l+i?“2)e''o^4(^+“h(2-a)(';^_^^gl2.)-2^)^^_^ |4't|^+ max |VT|2+ max IV^T]), 

B+x[0,r] S+x[0,r] B+x[0,T] 

(5.22) 


with C > 0 independent of i?, T, 0. Therefore, 


Cw < MgiL(|Vp|). 


(5.23) 


(c) Let w = w for A,/r > 0. Then 

Cw = Cw + \C{e-p^-) < M^k{\Vp\) - Xp^Anne-P^-^ - Xpke-P^ 
< MgK(|Vp|) - Xp{c2P - c^)k{\Vp\)e-P^-. 






















Global Estimates for Generalized Forchheimer Flows 


19 


Choose /i = (1 + C 4 )/c 2 and 


A = 


(5.24) 


Then on x (0,T], Gw < 0. Since p{x,t) G C^’^{U x [0,T]), the function K{\Vp\) is bounded 
below by a positive number. Then by (j5.6j) . £ is a parabolic operator. Therefore, the maximum 
principle for operator C implies 


max w= max w. 

B+x[0,r] 9p(B+x(0,T)) 

Here dp denotes the parabolic boundary. 

When (x,t) G {B^ x {0}) U {{dB^ \ Tij) x (0,T]), we have Cix,t) = 0 hence 

w{x, t) = 0 + < A. 

When (x,t) G T^j x (0,T]: v{x,t) = 0 hence w{x,t) = A. 

Thus, 

max w = max w<X = w(x,t), V(x,t) G Tr x (0,Tl. 

B+x[0,T] dp{B+x{0,T)) 

Hence, we have Wx„ < 0 on T/j x (0, T], equivalently, Wx„ < pX on T^j x (0,T]. Note that w = v on 
r_R /2 X [0T,T]. Thus, on rj:j /2 x [6T,T] we have 

KoiPxr. -^xj = Vx^ = Wx^ < pX, 

which implies 

Px^ < + Kq VA. (5.25) 

Replacing p by —p, and by — we obtain again (I5.25P with —px^ —^xn place of Px^ and 
* T'herefore^ ^ve have on rfl /2 X [6»r,r] that 

\PxJ < C{\^xJ+l^o^pX). 

Combining with the fact pxi = ^xi on Tr x (0, T] for i < n, we assert that 

|Vp| < C(|VT| + Kq VA) on rjj /2 X [eT,T]. (5.26) 

Using (15.241) and (I5.22p in (I5.26p . we obtain 

max |Vp(x, t)| < (7(1 + '^(1 + (0r)“^)^° exp (C' max |p|) 

rfl /2 X [eT,T] B+ X [eT/2,T] 

2 

•(l+_max IT^I 2^ + _ max |VT|^+_max iV^Tj), (5.27) 
B+x[0,T] B+x[0,T] B+x[0,T] 

thus proving (15.9p . □ 

The general domain and boundary are treated in the next theorem. 

Theorem 5.2. Let p{x, t) be a solution of (j3.ip . Then for any Tq > 0, T > 0, and 9 G (0,1), 


max |Vp(x,t)l < (7(1 + (0T) ^)^°exp(C'_ max Ipl) 
rx[To+6»r,To+T] Ux[To+eT/2,To+T] 

2 

• (1 + _ max + _ max |V4'p + _ max |V^4'|), (5.28) 

Ux[To,To+T] Ux[To,To+T] UxITo,To+T] 


where C, C are positive numbers independent ofT^, T, and 9. 
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Proof. By replacing p{x, t) with p{x, Tq +1), we can assume, without loss of generality, that Tq = 0 
and p € x [0,T]). Let xq G dU. There exist an open neighborhood V of xq, a radius i? > 0 

and C^-bijections y = <l>(x) : U CiV ^ and x = T(y) : U (IV such that = T“^, 

<l>{TnV) = B2Rn{yn = 0}, (5.29) 

and 

ll^llc'2(c/nv')> ll^llc2(B+j) - ^0 <^0 > 0. (5.30) 

Define p{y) = p(T{y)) and 'l'(y) = T(T(y)). We use the fact that p{x,t) is a solution of (15.ip . 
Simple calculations give 


Pxi — ^^Pyk^Xii PxiXj — PvkVl^Xi^Xj + 'y^.Pyk^xja 

k,l=l 


k=l 


k=l 


with = {^^)k=i,...,n- Thus, 


Pt-Y^ ^ijPx^Xj =Pt-Y -^klPym + Y ~^^Pyk = ^P' 

i,j=l k,l=l 


k 

XiXj • 


fc=l 

where 

n n 

Aki = Y ~^k = -Y 

i,j=l i,j=^ 

Therefore, p{y, t) satisfies the equation 

£p{y,t) = 0 on X (0,r), p{y,t) = ^{y,t) for = 0. 

Now we check the conditions (j5.5h . (15.6h and (15.7h . By chain rule and (I5.30h . 

|Vyp(y,t)| < co\Vp{x,t)\, \Vp{x,t)\ < co\Vyp{y,t)\, 

IV^TI < codV^Tl + |VT|), IV^TI < codV^Tj + \Vy^), 

for some cq > 1. We have for any k,l = 1,2,... ,n that 

n 

\^ki\ = I Y ^ cl max \Aij\ < n^cl{l + a)d2{l + |Vp|)“ 

*j=i 

< n^Co(l + a)d2{l + Cq dVj^p|)““ < n^Co(l + a)d2Co(l + \Vyp\)~°'. 
This yields ()5.5p . For all ^ G M"’, 

n n n n 


Y^ ■^klf.kf.l — Y^ Y1 ■^d^Xi^Xj^kf.l — Y^ AijrjiPj 

k,l=li,j=l i,j=^ 


k,l=l 


with rj = (D<h)^, i.e., ^ = {DT)p. By (15.3h we have 


(5.31) 


(5.32) 


(5.33) 


(1 - a)di{l + |Vp|) < Y •^ij'niPj < d 2 {l + |Vp|) 

*d=i 


— CL I 2 


(5.34) 
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Note that Cq < |??| < co|^|. It follows from ()5.34p that 


•AkiCkCi<cld2{l + Co^\Vyp\) “ICP < cgd2c;5(l + |Vyp|) “lip, 

k,l=l 


^ AkiCk^i > (1 - o)co ^di(l + co|Vyp|) “lip > (1 - a)cQ “(1 + \Vyp\) “||p, 
k,l=l 

hence ()5.6p is satished. Next, we bound | 6 fc| for each /c = 1, 2,... , n by 

n 

\bk\ < \ < con2(l + a)d 2 (l + |Vp|)"“ 

*j=i 

< con^(l + a)d2il + Cg pVyp|)““ < con^(l + a)d2Co(l + |Vyp|)““. 

Hence we have p5.7p . 

Applying estimate p5.9p from Proposition 15.11 to p, ^ and equation p5.3ip . we obtain 

max \Vyp{y,t)\ < Cr{ 1 + {9T)~^)^° ex.p (C' max |p —^|) 
rfl/ 2 x[ 6 »T,r] B+x[eT/2,T] 

•(l+_max |^tP°+_max |Vy^p+_max |V^^|). (5.35) 

B+x[0,T] B+x[0,T] B+x[0,T] 

Let p > 0 sufficiently small such that U H Bp{xo) C Thanks to ^5.291) . (I5.35P and ^5.321) . 

it follows that 

max |Vp(x, t)| < Ci?(l + exp (C'_ max |p|) 

(rnBp(xo))x[er,T] ux[eT/2,T] 

■ (1 + _max |'I'fP° + _max iV'kp + _max (IV^TI + IVTI)). (5.36) 

Ux[0,T] ;7x[0,r] i7x[0,T] 

By Cauchy’s inequality for the last |VT|, we obtain 

max |Vp(x, t)| < Cr( 1 + (0r)“^)^° exp (C'_ max |p|) 

(rn.Bp(a;o))x[er,T] Ux[eT/2,T] 

■ (1 + _max iTiP^ + _max iVTp + _max IV^TI). (5.37) 
(7x[0,r] Ux[0,T] ;7x[0,T] 

By using a finite open covering of T, we obtain the desired estimate p5.28l) from p5.37p . The proof 
is complete. □ 

The bounds of ||Vp(t)||ioo(r), in fact, can be expressed in terms of the initial and boundary 
data as follows. 

Corollary 5.3. Let p{x,t) be a solution of (13.ip . and let a satisfy (12.1511 . 

(i) If 0 < t < 3 then 

||Vp(t)||ioo(r) < Ct-^<^il + max(||TP|^°^ + ||VT||ioc + HV^vkHi^)) 

[t/4,Ej 

•exp|c't"^(l+||po||L“(i7))"''^(l + [^nuA(Q;,t)]“)'”^(l + ||VT||2,<^9i(;7x(0,t)) + II^PlL“U(C/x( 0 ,t)))''^}- 

(5.38) 
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If t > 1 then 

||Vp(i)||ioo(r) < C(1 + + ||V^||ioo + HV^M-Hi^)) 

[t l,tj 

■exp ^C'{l+\\po\\La(^^^^y^{l+[EnvA{a, i)] “““)^®(l+||V'I'||ia9i([/x(t-i,t)) + ||'I't||i,“9i((7x(t-i,t)))^^|' 

(5.39) 


(ii) If A{a) < oo then 

limsup ||Vp(t)|| 2 ,oo(r) < C{1 + limsup(||^'i||^“oo + ||V^'||ioo + ||V^^'||loo)) 

t—>-oo t—>-oo 

•exp|c''(l + ^(a)“)^3(l + limsup(||V^'(t)||i«9i(;7x(t-i,t)) + ll^i(i)llL“«(c/x(t-i,t))))''"|- (5-40) 

t^oo > 

(iii) If I3{a) < oo then there is T > 0 such that for all t >T 


l|Vp(t)||Loo(r) < C(1 + maxdl^tllf^ + ||V4/||i^ + HV^M-ld^)) 

[t l,tj 


•exp 




(5.41) 


Proof, (i) For 0 < t < 3, applying Theorem 15.21 to Tq = t/4, T = 3t/4, 6 = 1/3 and using estimate 
(I4.22P we have 


l|Vp(t)||Loo(r) < Ct->^yi + max(||Td^°c« + ||VT||ioo + HV^THloo)) 

[t/4,i] 

• exp ^C't (1 + ||po||l“([/) + [EnvA{a, t)] a-a)^3(l + ||VT|dQ9i((7x(o,t)) + l|4'4||L“'Ji(C/x(o,t)))^^|' 
Then (|5.38p follows. 

When t > 1, applying Theorem 15.21 to Tq = t — 1, T = 1, 0 = 1/2 gives 
||Vp(t)||ioo(r) < (7(1 + max (||Tt||^^ + ||VT|||oo + ||V^T|doo)) • exp |c' max ||p||l°°|. (5.42) 

Using estimate (j4.23p we have 

||Vp(t)|doo(p) < C(1 + max(||Tt||^°^ + ||VT||i^ + HV^TlUco)) 

[4-1,t] 

• exp |c''(l + ||po||l“((7 ) + [EnvA{a,t)] “-“)^®(1 + ||V'I'||£,a9i(;7x(4-i,t)) + ||'l't||i,«9i([/x(4-i,t)))^^|- 
Then ()5.39p follows. 

(ii) If 74(a) < oo then taking limit superior of (|5.42l) and using ()4.24l) yield 


limsup ||Vp(t)||poo(r) < C{1 + limsup(||Tdp?o + ||VT|||oo + ||V^T||l<x>)) 

t^OO t^oo 

■ exp |C''(1 + A{a)°‘-<^y^{l + limsup(||V'I'(t)||£,t79i(f/x(t-i,4)) + l|4'4(^)llL“u(t/x(4-i,t))))'^^|- 


Thus we obtain (15.400 . 
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(iii) Using ()4.25l) in ()5.42l) gives 


||Vp(i)||ioo(p) < C{1 + rnaxmtWiL + ||V4/||ioo + 

[t l,t\ 


•exp 


for alH > T with some T > 0. This proves (I5.41h . 


VT| 


L°“ii{Ux{t-l,t)) + ll^i||L“5i(C/x(t-l 




□ 


6 L®-estimates for the gradient 

In this section, we estimate the pressure gradient in L^-norm for any 0 < s < oo. Throughout 
this section, p{x,t) is a solution of (13.ip . First, we establish the basic step for the Ladyzhenskaya- 
Uraltseva iteration. 

Lemma 6.1. Let s > 0, Tq > 0, and T > T' > 0. Define 

Mf, = max iVpl and u = max{| VpP — M?, 0}. 
rx[To+T',ro+T] 

Let C(t) be a smooth cut-off function on [ro,To + T] with (^ = 0 on [TojTo + T']. Then 


sup / v^~^^{x,t)Cdx-\- 
[ 7 o,To+r] Ju 


'To Ju 


K{\\Jp\)\\J^p\^v%dxdt < C 


I To Ju 


v^~^^\C,t\dxdt. (6.1) 


Proof. Without loss of generality, assume Tq = 0. Note that uC = 0 on T x [0,T]. Denote 
X = X{i;>o}- Multiplying the equation (13.11) by —V • {v^(Vp), integrating the resultant over U and 
using integration by parts, we obtain 


1 d r 

2s T 2 dt Jjj 


C,dx = — I dj{K{\Vp\)dip) di{v^djp()dx-\- 
i,j=^ 


1 


lu 


2s+ 2 


[ v''~^^Ctdx. (6.2) 

Ju 


Using the product rule for the first term in the right hand side of (j6.2l) (see Lemma 3.6 of m) we 
rewrite above equation as 


1 d 


25-1-2 dt Jjj 


L 




dyi{K{\y\)yi 


n n 


-»E E 

i,j=l l,m=l ' 


U 


y=Vp 

s-1; 


djdip djdipv^Qdx 


dyi{K{\y\)yi) djdip djp{v^ ^didmpdmPx) Cdx+ \ [ 

y=Vp J Is I Jjj 


v^+^Ctdx. 


u 


We denote the three terms on the right-hand side by Ii, I 2 , and Is. It follows from the calculations 
in Lemma 3.6 of m that 


h<-il-a)Y, / i^(|Vp|)|V(9,p)|VCdx, 


j = ,JU 


/2<-(l-a)s^iF(|Vp|)|v(i|Vp|- 


'U 

Combining these estimates, we find that 
1 d 


^xCdx < 0. 


+ (1 — a) [ K(\Vp\)\V‘^p\‘^v^Cdx <—^— [ v^^^Ctdx. (6.3) 
Ju 2s+ 2 Ju 


2s T 2 dt Ju 

Inequality (|6.1I) follows directly by integrating (|6.3I) from 0 to T. 


□ 
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In the following proposition, we iterate the inequality in Lemma l6.II in order to estimate W^’^- 
norm of p in term of its and L°° norms. 

Lemma 6.2. Let s > 1, Tq > 0, T > 0, and 0 < 6' < 6 < 1. Define 


Mfe = max llVpILoocp') and v = max{|Vpp — M?,0}. 
[To+e'T,To+T] ' 


Then 


rTo+T 


s—1 pTq-\-T 


j K{\Vp\)v^dxdt<CTU{s) + Cdl-’^ j j {l + \Vp\^-^)dxdL (6.4) 
Ju JTo+e'T Ju 


JTo+eT JU 

where constant C > 0 is independent oJTq, T, 9, and O', 


U{s) 


0, «/ 1 < s < 3 — o, 

s-l I 1 

do + ^, if s > 3 — a, 


(6.5) 


do = NQ{{e - e')T) ^ + with No = sup ||p||l°°({/)- (6.6) 

[To+e'T,To+T] 

Proof. Without loss of generality, assume Tq = 0. The proof consists of three steps. 

Step 1. Let C{t) be the cut-off function with ^ = 0 for t < O'T and () = 1 for t > 6T. For 
•s > 0, by applying Lemma [2]3] with k = 0, and s -|- 1 in place of s, multiplying (12.231) by C^(^) and 
integrating from 0 to T, we have 


'0 Ju 


K{\Vp\)v'^^‘^C^dxdt < Cmaxlpl"^ f f K{\Vp\)\V‘^p\‘^v'^C^dxdt + CM^ f f K{\Vp\)v'^(‘^dxdt. 

Jo Ju Jo Ju 


To estimate the first integral on the right-hand side, we use (|6.ip with in place of C, and find 
that 



K{\Vp\)v^+^C‘^dxdt<CN^ [ [ v^+\Ctdxdt + CM^ [ [ 

Jo Ju Jo Ju 


K{\Vp\)v^C^dxdt. 


(6.7) 


Now using Young’s inequality we have 


K{\Vp\)v^+^C^dxdt<- I I K{\Vp\)v^+'^Cdxdt + CN^ I I K-WVp\)v^Ctdxdt 


s-\-2/~2 


10 Ju 


10 Ju 


>0 Ju 


+ 


CM^ f [ K{\Vp\)v^C^dxdt. 
Jot Ju 


Thus 

rT 


[ [ K{\Vp\)v''+^C^dxdt<CN^ 

Jo Ju 


10 JU 
Note that 


K 


>0 Ju 


-^{\Vp\)v''Qdxdt + CMl [ [ K{\Vp\)v%^dxdt. 

Jot Ju 

( 6 . 8 ) 


K-\\Vp\)v^ < CiL(|Vp|)(l + \Vp\f%^ < CiL(|Vp|)(l + {v + 

< C7L(|Vp|)(l + < CiL(|Vp|)(u"+“ + 1 + 

and 

K{\Vp\)v^ < CiL(|Vp|)(l + u"+“). 
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Thus (j6.8l) implies 

' K{\Vp\)v^+\^dxdt<C[N^{{9-e')T)-‘^ + M^] Cl K{\Vp\)v^+-X{c>o}dxdt 

Jo Ju 

+ ' 


10 Ju 


c r [ [N^{{e-e')T)-\l + M^^^+^^) + M^]xic>0}dxdt. 

Jo Ju 


Therefore, 


r/ 

lex Ju 


K{\Vp\)v^+‘^dxdt<C[N^{{e-e')T)-"^+ M^] ! [ K{\Vp\)v^+^dxdt 

Je'T Ju 

+ CT[N^{{9 - e')T)-^{l + + M^]. 


Let Ni = TN^{{e - 6l')T)-2 and N 2 = Tdo- Then 


r/ 

I ex Ju 


K{\Vp\)v^^^dxdt <Cdo [ [ K{\Vp\)v^+'^dxdt + + CN 2 . (6.9) 

Je'x Ju 


Equivalently, for s > 2 

rX 


lex Ju 


K{\\/p\)v^dxdt <Cdo [ [ iL(|Vp|)u"-(2-“)dxdt + CiVi(l + Mfc2)^-(2-a) ^ 

Je'x Ju 


Step 2. We prove (16.41) for s G [1, 3 — a]. First, 


( 6 . 10 ) 


fl 


K{\Vp\)vdxdt< [ [ K{\Vp\)\Vp\‘^dxdt < [ [ {1 + \Vp\‘^-^)dxdt. (6.11) 

Jex Ju Jex Ju 


This yields (|6.4p for s = 1. 

Second, let s = 1 — a in ()6.8I) . 


< 


[ K{\Vp\)v^-'^dxdt < CN^{{9 - e')T)-^ [ [ (1 + |Vp |)[ [ v^-'^dxdt 

Ju Je'x Ju Je'x Ju 

CN^{{9 - 9')T)-^ [ f {l + \Vp\y\Vp\^^^-^Uxdt + CM^ f [ \Vp\‘^^^-^Uxdt. 

Jo Ju Jo Ju 


lex JU 


Therefore, 


[ [ K{\Vp\)v^-^dxdt<Cdo [ / (1 + |Vpp-“)dxdt 

Jex Ju Jo Ju 


( 6 . 12 ) 


This implies (16.41) when s = 3 — a. 

Third, when s G (1, 3 — a), let fJ be the number in (0,1) such that ^ + 3^. Then, by 

interpolation inequality: 


ex Ju 


K {\V p\)v^ dxdt] < 


ex JU 


1-/3 


K{\Vp\)v'^-‘^dxdt 


ex JU 


Using ()6.12p to estimate the last double integral, we obtain 


11 

lex Ju 


K{\Vp\)v^dxdt < (^j j {I + |Vpp-“)dxdt) 


s{l-J+yh) 
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Note that j5s/{‘i — a) = (s — l)/(2 — a). Then 

rT r 


s-i rT 


IdT JU 


K{\Vp\)v^dxdt < Cdg 


10 Ju 


(1 + \Vp\^-'^)dxdt. 


(6.13) 


Therefore, (16.lip and (I6.13P imply (16.41) for s E [1, 3 — a]. 

Step 3. When s > 3 — a, let m € N such that 

s — (3 — a) s — (3 — a) 

< m < —^- - + 1, 


2 — a 


2 — a 


then s—m{2—a) E (1, 3—a]. For each k = 1,2,... 6^ = 9—{d—d')k/m and Qk = Ux {OkT, T). 

Also, let Ckit) be a smooth cut-off function which is equal to one on Q^. and zero on Qt \ Qk-i- 
There is a positive constant c > 0 depending on U, such that < mc[{0 — 9')T]~^ , for all 

A: = 1, 2,... , m. 

From (j6.inp we have 

T T 

[ [ K{\Vp\)v^dxdt < Cd^ [ [ K{\Vp\y-^^^-°^'>dxdt 

Jot Ju Jb^t Ju 

+ C*A^2(1 + do T <^0 T ... + dj^ ^). 

Using inequality z* < (1 + z)'^~^,\/z > 0, we infer 

T T 

[ [ K{\Vp\ydxdt < CdJJJ [ [ K{\Vp\y-^^‘^-°-^dxdt 

Jbt Ju JomT Ju 


+ [1 + + CTdo{l + do)™“^- 


Therefore we obtain 

rT 


[ [ K{\Vp\)v^dxdt < Cd'^ [ [ K{\Vp\)v^-^^^-'^Uxdt 

JdT Ju JdmT Ju 

+ + CTdoil + do)”'"^ 

Since s — m{2 — a) E (1 — a,3 — a], estimating the double integral on the right-hand side by 
p6.13p gives 


\m—1 


[ [ iF(|Vp|)u"dxdt <CiViM^^^"“^^"“^^[l + doM“^^^““^]™-i + Crdo(l + do)^ 

Jot Ju 

s-m( 2 -a)-l rT r 

TCdg / / (1 + |Vp|2-“)dxdt. 

Je-mT Ju 

Since m < and then 

Z — CL 2i — CL L — (1 ' 

[ [ K{\Vp\)v^dxdt < + Crdo(l + do)^"^ 

JdT Ju 

+ Cdf=^ [ /(IT |Vp|2-“)dxdt < + d^~^M^) + CT{do + d|^) 

JOmT Ju 

+ Cdf^ [ [ (1 + |Vp|2-“)dxdt. 

JOmT Ju 
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Note that 


iVi(M,+ 4 -“ ^Mi) + T{do + ) < r(dodf + dod^- + do + do^-“) 


-"0 
s-l I 1 


= r(v“ + ^2-+^ + dp + / 2 -a^ < cr(v“ + ^2-^+i + < CT(do + '"), 


since 1 < ^ for s > 3 — a. So, 

cT 


[ [ K{\Vp\)v^dxdt < CTU{s) + Cd^^ [ / (1 + \Vp\‘^-‘^)dxdt. 

JeT Ju Je'T Ju 

This completes the proof of (|6.4p for all s > 1. 

From the estimates of K(\Vp\)v^ is Lemma 16.21 we derive direct estimates for |Vp|‘ 

Proposition 6.3. Let Tq > 0, T > 0 and 0 G (0,1). If s > 2 — a then 


rTo+T 

Jto+OT 

and if s > 2 then 


[ ivpr 

Ju 


dxdt < CT{1 + T>") + CV^ 


--1 


rTo+T 


/To Ju 


|Vp| °‘dxdt, 


sup 

[To+eT,To+T] Ju 


where eonstant C > 0 is independent o/Tq, T, and 9, 


— — 1 


fTo+T 


\Vp\^dx < ce-\l + V^) + C{0T)-^V^ 


/To Ju 


(6.14) 

□ 


(6.15) 


|Vpf ‘^dxdt, (6.16) 


max{^^, 2 ^ “!}) if‘I — a<s< 3(2 — 


s = 


S 

2-a’ 


if s > 3(2 — a), 


and 


V = {eT) 1 sup \\p\\loo^u)+ sup |Vp|2oc(r)- 

[To+eT/2,To+T] [To+eT/2,To+T] 


(6.17) 


(6.18) 


Proof. Without loss of generality, assume Tq = 0. Let Mb, v and do be defined as in Lemma [62J 
Proof of (|6.15p . From the definition of the function v and inequality \a — b\J > 2~JaL — for 
^ > 0,7 > 0 , we have 


/7 

ter Ju 


K{\Vp\)\Vp\^^dxdt< [ [ K{\Vp\)vMxdt + CTMiL 

JeT Ju 


We apply (|6.4I) with 0' = 0/2 to estimate the integral on the right-hand side. It results in 


[ [ K{\Vp\)\Vp\‘^^dxdt < CPU{s) + CTM^^ + Cdt'^ [ [ (1 + \Vp\‘^-'^)dxdt. (6.19) 

JeT Ju Jo Ju 

For s > 3 — a, then by definition ()6.5p and then Young’s inequality, it follows 

[ [ KilVpDlVpl'^Mxdt <CT{do + dl 

JeT Ju 


a-l , 1 S-l S-l rT 

^+df+ dl-'^) + Cdl’^ I 


|Vp|2-“dx 


s-l I 1 
2-0“'“ 2 ' 


/o JU 
-1 rT 


<CT{l + d^-'^ ')+Cd, 


/o Ju 


\\Jp\^-‘^dxdt. 
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Estimating i^dVpl) from below by (|2.11l) . we get 

rT 


[ + + [ [ \Vp\‘^-^dxdt. ( 6 . 20 ) 

Ju Jo Ju 


leT 


Replacing 2s — a by s > 2(3 — a) — a = 3(2 — a) in (16.201) we obtain 

CT s s — {2 — a) pT 


[ \Vp\^dxdt < CT{1 + [ [ \Vp\^-^dxdt. (6.21) 

JeT Jo Ju 


For 1 < s < 3 — a, we have from (|6.19l) that 

rT r 


s —1 

.s /2 I ^ 2 -a' 


JOT JU 

Using property ()2.11l) again yields 

rT 


K{\Vp\)\Vpy^dxdt < CT{d"^" + cio'-“) + Cd^ 


s-i rT 

72 — 0, 


/ / |Vp| 

Jo Ju 


2—a 


dxdt. 


[ [ \Vp\‘^^-^dxdt <CT{l + d"J^+ dy^)+ Cd, 
Jot Ju 

Replacing 2s — a by s G [2 — a, 3(2 — a)] gives 


«-i rT 
2 —a I 
0 




s + a s-{2-a) 

- 2 ( 2 - 0 ) 


\Vp\^dxdt < CT{1 + ) + Cd^ 


s—(2 — a) pT 
2(2-o) / 


/o Ju 


\Vp\^-^dxdt 


JeT Ju Jo JU 

<CT(l + do ^ ) + Cdo' ’ / / |Vp|2-“dxdt. (6.22) 

Jo Ju 

Note that < CD. Then combing (j6.21l) . for the case s > 3(2 — a), with (|6.22l) . for the case 
s G [2 — a, 3(2 — a)], we obtain (16.151) for all s > 2 — a. 

Proof of (I6.16h . For s > 2, we have from (16.Ih with O' = 6/2 that 


sup [ \\/p\Udx<C sup f {v^^^ + M§)Cdx<CM§+ [ f v^/^Ctdxdt 
[ 0 ,T] Ju [ 0 ,T] Ju Jo Ju 

<CMI + C r [ \Vp\^/^\Ct\dxdt. 

Jo Ju 


Choose Q{t) such that (" = 0 for t < 3dT/4, and ^ = 1 for t > 6T. Then 

sup [ \Vp\%dx + C{OT)-^ [ [ |Vp|*dxdt, 

WT.T]Ju J^tJu 


(6.23) 


with do being defined by using, again, 9' = 9/2. We apply (16.151) with 9 being replaced by di = 3d/4 
and do, hence V, dehned by using 9' = 29\/‘J>. We obtain 


3-1 


< C9-\l + V^) + C{9T)-^D-- 


sup [ \Vp\^dx <CD^/^ +09-^1+ D~^) + C{9T)-^DT^-^ [ [ \Vp\ 

[eT,T] Ju Jo Ju 

rT r 

\ —I'T'i.TrT' —1 / / IV7»,|2—a. 


2—a 


dxdt 


>0 Ju 


iVpU'^dxdt. 


Thus, we obtain (j6.16j) . This completes the proof of the theorem. 


□ 
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Now, we combine Proposition 16.31 with estimates in section [3] to express the bounds for the 
gradient’s L^-norms in terms of the initial and boundary data. 

In the following a always satisfies ()2.15l) . Let pi, qi be fixed and (5i be defined as in Theorem 
HU Also, let zi, Z 2 , Z 3 be defined as in Corollary 14.31 and /tq be defined by (j5.inp . We define a 
number of constants and quantities that will be used for the rest of the paper; 

Ki = max{2/r0) 1 + 2/5i}, (6.24) 

/Cl (t) = 1 + [EnvA{a, t)] , JCi{t) = 1 + /3{a) “-2“ + sup A(a, (6.25) 


/C2(t) = 1 + sup(||T||loo + + ||VT||ioo + IIV^TIIlco) 

[04] 


+ (I|V'I'IIl“«i(;7x(04)) + II'^'4|U“«(c/x(04)))^^) (6.26) 


^ 2 {t) — IT sup (||'I'||l°° T T llVTlI^oo T llV^TlIx,^^) 

[t-2,t] 


+ (l|V'^'llL“9i([/x(t-2,t)) + \\'^t\\L°‘ii{Ux{t-2,t))y^ ■ (6.27) 


Theorem 6.4. Let s > 2. 

(i) If 0 < t < 3 then 


'U 


\Vp{x,t)\^dx < T ||po|U«(t7))2-3+2/Ci(t)2-3;C2(t)2^'(l T rGi(T)dT) 


• exp 


|cT + ||po||L°i(f/))^®/Ci(t)^®(l T ||V'I'||i^a9i(i7x(o,t)) T ||'I't||i,“9i([/x(o,t)))^^|- (6.28) 


(ii) If t > 2 then 


'U 


\Vp{x,t)\^dx < C(1 T ||po||L<^(c/))'*"^+“/Ci(t)2-3+«^2(t)'"(l T £^Gi(r)dr 


• exp ■ 


|C''(1 T ||po|liQ([7))^i(i)^®(l T ||V'I'||i;^Q,9i((7x(i_2,t)) T ||'I't||i,“9i((7x(t-2,t)))^^|- (6.29) 


Proof, (i) Let t E (0, 3]. Applying (|6.16p to Tq = 0, T = t and 0 = 1/2 gives 


’U 


\Vp{x,t)\^dx < C(l + V~^ + Ct-^V^-^ 
where V is defined by (16.181) . Since s > 2 ^ — 1, then 


\Vpf-^dxdT), 


0 JU 


lu 


\vp{x,t)\^dx < ct-yi + vy{i + 



0 Ju 


|Vp| °‘dxdTj, 


(6.30) 


V<Ct ^ sup IIpIIIoc T sup ||Vp||^oc(r) 

[i/ 44 ] [i/44] 

< Ct~^ sup WpWl^ +Ct~^ sup ll^illloo T sup ||Vp|||oo/r). 

[i/ 44 ] [i/ 44 ] [i/ 44 ] 


Note that 
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To estimate ||p||l°o we use (|4.22l) . To estimate ||Vp||^oo(r), we apply (|5.38p of Corollary 15.31 We 
obtain 

V < + t-^Rl + = CV, 

where 

-^1 = (1 + l|Po||L“)^^^l(i)^®(l + l|V'I'||iQi9i((7x(0,t)) + ll'^'t|lL“5i(;7x(0,0))^^’ 

i?2 = sup ||4'||i,<x>, ^3 = 1 + sup(||Tt||^°oo + ||VT|||,oo + ||V^T||loc). 

[o,t] [o,t] 

By relation (I2.13h . we have |Vpp““ < <^(1 + H{\Vp\)), hence we can use (13.8h to estimate the last 
integral in ()6.30p . Substituting these estimates into P6.30D gives 


|Vp(x,t)|*(ix < Ct ^VHl + \\po\\l2+ I Gi{T)dT 


'U 


2 

L2 


f 


(6.31) 


hence 

[ \Vp{x,t)\^dx <Ct-^-^^^{Ri + R 2 + R 3 f^R 4 e^'^~^^''^\ (6.32) 

Ju 

where R 4 = (1 + ||po||l2)^(1 + fo Gi(T)dT). For the sum Ri + R 2 + R 3 , we clearly see that 
R 2 ,R 3 < K^ 2 {t), and 

Ri<c{i + \\po\\L<^^u)ric3{tric2it). 

Therefore 


\Vpix,t)\^dx < + IIpo||l 2)'(1 + ||po||L-)''"^[/Cl(t)^^/C2(t)]2^'(l + I^Gi{T)dT) 

• exp (6.33) 

Then by Holder’s inequality, we can combine the powers of and L" norms of po, and obtain 

(iraii . 

(ii) Let t > 2. Applying (I6.16P with Tq = t — 1, T = 1, 0 = 1/2, 


lU 


\Vp{x,t)\^dx < C{1 + V'^) + CV-- 


—1 


it-i Ju 
< 


\Vp\^-'^dxdT 


C(l+P)"[l+ f [ |Vp|2-“dxdrj, (6. 
L Jt-i Ju 


34) 


where P is defined in (I6.18p . We have 


12 


P < C sup IIpIIloc + ||Vp||£,oo(r) 


< C sup ||p||l°° + ||^||l°° + ||Vp||Loo(r) 


(6.35) 


Similar to part (i), but using (|4.23p and p5.39p instead of (|4.22p and (|5.38p . we have 

P < C{Rl + Rl + R^e^'^^), 


(6.36) 
































Global Estimates for Generalized Forchheimer Flows 


31 


where 


-^5 — (1 + \\Po\\l<^Y^ + \\'^'^\\L°‘n{Ux{t-2,t)) + \\'^t\\L°‘ii{Ux{t-2,t))Y^ ^ 

i?6 = sup 11'h 11 £,00 , Rj = 1 + max (H'l'i + || V'J'H^oo + || V^\h||£oo). 
Combining (|6.34p . ()6.36p and p3.10p gives 


J \\/p{x,t)\^dx < C{R\ + -Rg + R^e ^'+ llPollia + EnvA{a,t)°‘-<^ + Gi(r)dT^ 

<C{R5 + Re + R7e^'^Y^H^ + \\Po\\LY"}Ci{tr{^+ t Gi(r)dr). (6.37) 
Note that Rg,R 7 < Y 2 Y) and 

RG<c{i + \\Po\\LY"^iCi{tY^Y2{tY^. 

Therefore, similar to part (i), we obtain (|6.29l) from (I6.37p . □ 

Remark 6.5. (a) Regarding the regularity requirement for the initial data po{x) in estimates of 
||Vp||£s in later time t > 0, the right-hand side of (j6.28p and (j6.29p only need ||po||l“- This shows 
the regularity gain of the solution when t > 0. 

(b) The -estimates established above hold for s > n, hence the Holder estimates follow 
thanks to Morrey’s embedding. This is opposite to Ladyzhenskaya-Uraltseva’s proofs in im which 
establish Holder continuity first and use it to obtain -estimates, but not for all s. The reason 
for this simplification is our equation’s special structure, see (I2.8p ~ (l2.13p . 

For large time estimates, the bounds in Theorem 16.41 can be established to be independent of 
the initial data as in the following. 

Theorem 6.6. (i) If A(a) < 00 then 


limsup / \Vp{x,t)\’^dx < 

i—>-00 Ju 

• exp + limsup(||VT(t)||£a5i(f;x(t-i,i)) + ll^i(0llL“5i(c/x(i-i,t))))'^^|) (6.38) 


where 


.4,1 = 1 + ^(a) , ^i = l + limsup f Gi(r)(ir, 

t^oo Jt—l 


(6.39) 


.42 = 1 + limsup(||T(t)||£cx> + W^tiYwi: + ||VT(t)||ioo + ||V2T(t)||£^) 

t—>-oo 

+ limsup(||VT(t)||£Q9i([;) + tiYW 1°“^-{u)Y^ ■ (6.40) 

t—^oo 

(ii) If I3{a) < 00 then there is T > 0 such that for all t > T, 


lu 


\Vp{x,t)\’^dx < CICi{tY'^^’>+^]C 2 {tY" (l + ^ Gi{T)d7 

|c''^l(t)^®(l + ||V'I'(t)||£a9i((7x(i_2,t)) + {u x(t-2,t)}Y‘^"^ ■ (6-41) 


• exp 
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Proof, (i) Taking the limit superior as t —)• oo in (|6.34l) and (I6.35P gives 

limsup / \Vp{x,t)\^dxdt 
t^oo Ju 

< C'(l + limsup(||p||ioo(f/) + ||T||ioo(^) + ||Vp||ioc(r)))^^ • (1 +limsup [ f \Vpf~^dxdT). 
t—>co t—>-oo Jt—1 Ju 

(6.42) 

Applying (|4.24l) . (|5.40l) and (I3.12p gives 


limsup / IVp(x, 
t—>-oo Ju 


< C 


+ limsup(||V'I'(t)||ia9i (ux{t-i,t)) + ll'I't(^)llL“9i((7x(t-i,t)))) 


22 


t—>-oo 


2 

+ limsup 114'11 LOO (ij) + limsup(l + + HV'kH^oo + || V^THloo) 


t—>-oo 


t^oo 


exp{C'.Ai®(l +limsup ||4't(t)||L-w({/x(t-i,t)))^^} x (l + ^? +limsup [ Gi(r)dr). 

t^OO -I X t^oo Jt—1 ' 


Then estimate (j6.38h follows. 

(ii) By combining (|6.34p . (|6.35p with (14.2511 . (I5.4ip and (j3.14p . we have 


/ \Vp{x,tWdx 

Ju 

< c|[l + /3(a)“-2“ + ||A(a, ] ^[1 + \\^^\\L°‘nux{t- 2 ,t)) + \\'^t\\L°‘<iiux{t-2,t))Y 

L L“-“ {t—2,t) \ V w/ 

+ sup 11T11 Loo (f/) + (1 + sup + llVTlI^oo + IIV^TIIloo)) 

■ expIC')! + /3(a)“-2“ + \\A{a, •)|| “IS_ ] ^[1 + ||V 4 '||L 09 if/x(t- 2 ,t)) + \\^t\\L°‘nux{t -2 

La-a (t—2,t)-' \ ’ n \ 


, 0 ) 1 '"}} 

(l + /?(a')“-2“ + /1 (q:, i — 1 )“-“+^ Gi(r)tfr) 


for alH > T with some T > 2. Note that \\A{a. dH ° a < 2^/“ supfj _2 a A(a, Then 

L°‘-a(t-2,t) ‘ 


|Vp(x,t)|®dx 

+^2(t)exp{C"^i(t)^3[i + 


+ ll^t||L“9ii7x(t-2,t))]^^}| 


Therefore, inequality (I6.4ip follows. □ 

Remark 6.7. In case T = T(x), we want to investigate the long-time dynamics by using the notion 
of global attractors, see, e.g., However, the independence of estimates (I6.38P and (I6.4ip 

on the initial data does not yet prove the existence of a absorbing ball. Therefore the existence of 
the global attractor is still an open problem. 
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7 L°°-estimates for the gradient 


In this section, we obtain global L°°-estimates for the gradient of pressure. Let p{x,t) be a solution 
of IBVP (13.ip . For each m = 1, 2,..., n, denote Um = Pxm u = {ui,U 2 , ■ ■ ■, Un) = Vp. We have 


^ = dm{y • {K{\u\)u)) = V • {K{\u\)dn,u) + V • 


(7.1) 


Since diUm = dmUi, we have dmU = (dmUi, , dmUn) = {dlUm, • • • , dnUm) = ^Um, and UidmUi = 
UidiUm = u ■ Vum- Therefore, we rewrite (17.111 as 


= V . (A'(|u|)Vu„) + V . [k'{\u\)F^u 


We consider this as a linear equation for Um with variable coefficients, and write 

dUm 


dt 


= V • A{x,t,VUm), 


(7.2) 


(7.3) 


where 


A{x,tA) = K{\u{x,t)\)^ + K\\u{x,t)\)-^^-T-i—^uix,t)- 


\u\ 


Using properties (I2.11|l and (I2.12|l . one can verify that 

|^(a;,f,OI < (1 + a)K{\u{x,t)\M\, 
A{x,tA) A > (1 - a)7sr(|u(x,f)|)|^|^ 


(7.4) 

(7.5) 


We will apply De Giorgi’s technique to equation (17.2p . In the following, we fix a number sq 
such that r = So satisfies (j2.19D . Note that Sq > 2. We will also use Sj for j > 1 to denote some 
exponents that depend on sq but are independent of a. Let 

si = (I-2/s*)-i > 1. (7.6) 

Theorem 7.1. For any Tq > 0, T > 0, and 9 G (0,1), ift£ [Tq + 6T^ Tq + T] then 

l|Vp(f)||Loo('[/) < G(1 + ( 0 r) ^) 2 A 2 ||Vp|h2([/x(To+0T/2,ro+T)) + C* sup ||Vp|hoo(r), 

[ro+6»r/2,To+T] 

(7.7) 

(7.8) 


where 


X = XiTo,T,e) = 


To+T 


(1 + |Vp|) dxdt) , 


lTo+eT/2 Ju 

and constant C > 0 is independent o/Tq, T, 9. 

Proof. Without loss of generality, assume Tq = 0. Denote Mi, = suprx[ 6 »T/ 2 ,T] I^pI- 
Fix m G {1, 2, ..., n}. We will show for t G [9T, T] that 

\\Px,n{t)\\L'^{u) < Cil + {9T) ^) 2 A 2 ||p3,^||^2(yx(To+0T/2,ro+T)) + 2Mb. (7-9) 

Let C(t) be a cut-off function with Q{t) = 0 for t < 9Tj^. We define for k > M^, u'm = 
max{um — 0}, then nXm = 0 on F. 
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Multiplying (I 7 . 3 jl by and integrating over U, using properties (j 7 . 4 h and ( 17 . 5 h . we obtain 

Integrating () 7 . 10 p from 0 to t for t G [ 0 ,T], and then taking supremum in t give 

max [ + C [ [ K{\u\)\Vu!;!^^\‘^Cdxdt < [ [ \ul^^\‘^CCtdx. ( 7 - 11 ) 

[o,T] Ju Jo Ju Jo Ju 

Let 

no = 4(1 - l/s 5 ) > 2 . ( 7 . 12 ) 

Applying Lemma 12.21 to function Um\ which vanishes on the boundary, weight W = K{\u\) and 
exponents r = so, g = g{so) = no, we have 




< c 


l^ikh^u) + {£ I^K{\u\)\V{ui^\fdxdty 
[ f Ar(|rt|) '^-‘odxdt 

*- ^0 J supp<^ 


ess sup \ \u^ 
4G[0,T] 

■T 


Z-SQ 


__£0_ gap 

Using () 2 . 1 ip . we have K{\u\) ^-“o < C{1 + |u|) 2 ““o, hence 


\u. 


L^oiQ,) < CAV-O 


max ; u 


(fc) 


-iG[0,T] 

By p 7 . 13 p . p 7 . 1 ip and the boundedness of function K{-) we find that 

rT 

lO JU 


C\‘^dx+ [ [ K{\u\)\V{ul^^C)fdxdtY. ( 7 . 13 ) 

Jo Ju ^ 


T 1 

L^o(Q,) < Xl“- 


( 7 . 14 ) 


Let Mo > Ml,. For i > 0, let ti = 6T{1 — ki = 2Mo(l — and let Q,, Q{t) be the same 

as in the proof of Theorem 14.11 For i,j > 0, define 

Ajj {(a;,t) . u,rni,x,t') > ki, t G (fj,T)}, Aj Aj j, 

Let Yi = ||l 2 (^. .). Applying p7.14p with k = fcj+i and C = Ci and then using the same 

arguments as in the proof of Theorem 14.11 (see detailed calculations in Theorem 5.6 of m) , we 
obtain 

Yi+i < for all i > 0, 

where n^ = I — 2/no, B = 4: and D = C(1 + 

We now determine Mq so that Lq < . This condition is met if 


Mo>C 


aV-o(i + ^)i/2j '/"Vo = CWi/2(l + 


Since Iq = \\um°'^\\L^{Ao,o) - \Wm\\L^{ux{eT/ 2 ,T)) and Mq > Mb, it suffices to choose Mo as 

Mo = CA^i/^(l + —) 2 \\Um\\L'2{Ux(eT/2,T)) + Mb- 


( 7 . 15 ) 
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Then Lemma [23] gives lirnTj = 0. Hence, dxdt = 0. Thus, Um{x,t) < Mq a.e. in 

U X {6T,T). Replace Um,u by —Um, —u and use the same argument we obtain 

\um{x.,t)\ < 2Mo a.e. in U x {6T,T). (7-16) 

By the choice of Mq we obtain from (17.1611 that 

\Umix,t)\ < C{1 + —) 2 X^^/‘^\\Um\\L‘^(Ux{eT/2,T)) + 

for all m = 1,... , n. Then (17.91) follows, and the proof is compete. □ 

We will combine Theorem O with the high integrability of Vp in section [6| to obtain the 
L°°-estimates. For the rest of this paper, we fix the following constants 

52 > max ( 2 , I and S 3 = 5i(2 - so)/so + 1, (7-17) 

t 2 — So J 

where sq and si are as in (|7.6I) . We use the same notation as in (|6.24p - ()6.27p . in Theorem 16.61 and 
also denote 

K2 = 1 + Si + S2K1S3, K 3 = S 2 Z 3 + a/ 2 , (7.18) 

^i(t) = 1 + /3(a)“-2“ + sup A(a, •)“^, (7.19) 


^ 2 {t) — 1 + sup (||T(t)||ioo + ||'I'l(t)||^(lo + ||V'I'(t)|||oo + ||V^'I'(t)||ioo) 

+ (l|V'^'llL“9i([/x(t-3,t)) + \\'^t\\L°‘ii(Ux{t-3,t))y^ ■ (7.20) 
We recall that the definition of s is given by (I6.17p . 

Theorem 7.2. (i) If 0 < t < 3 then 


IIVp(t) 1^00(11) < + 


\ S 3/2 


• exp 


Ct + \\'^'lt\\L°“ii{Ux(0,t)) + \\'^t\\L°“ii{Ux{0,t))y^ ■ (7.21) 


(ii) If t > 3 then 


||Vp(t)||ioc(c;) <C(l + ||po||L^(f/))"^^^/Ci(t)''3.3^2(i)^'^^^(l+ r Gi(r)dr)''^' 

t —2 

• exp 1(7(1 + ||po|l2a(i/))^l(0^^(l + \\'^'^\\L°“n{Ux{t- 2 ,t)) + \\'^t\\L°“ii{Ux{t- 2 ,t))y^'j ■ (7.22) 

Proof. First, we note in (|7.7p that 


O rTo+T I- ^-‘’0 

' (l + jVpiy^dxdt) ' 

TO+0T/2 Ju 7 


ll^7'llL2(f/x(To+er/2,To+T)) < ^ ^(1 + |Vp|)®^(ixdt 


To+eT/2 Ju 


1 

2 
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Hence, we have from (17.71) that 

sup ||Vp(t)|Uoo(^) < c(i + (0r)-i)^ ( [ 
[To+eT,To+T] '' Jt 


To+eT/2 Ju 


(1 + \Vp\y^dxdt 


*3/2 


+ C sup ||Vp||ioo(p). (7.23) 
[To+eT/2,To+T] 


(i) Let t G (0, 3], applying (j7.23l) with Tq = 0, T = t, and 6 = 1/2, we obtain 

l|Vp(t)||i,<x>([/) < [ [ {l + \Vp\y^dxdt] +(7 sup ||Vp||poc(r). (7.24) 

^ Jt/A Ju ' 


*3/2 


Using ()6.28p with s = S 2 and ()5.38p we obtain 

■ exp C't + \\po\\L°‘{i;)y^KLi{ty^{l + ||V'I'||pa9i([;x(0,t)) + ll'I't|lL“«i(!7x(0,d))^^ } 

+ Ct ^°/C 2 (/)exp C't +l|Po||l,“(i7))^®^l(i)^®(l+||V'I'||pagi([7x(o,t)) + ||^t||L“9i(;7x(0,d)) 

Hence, 

• {JC2{ty^^^+JC2{t))(l + Gi(r)dr)''^' 

C't + ||po||L“(i7))^®^l(i)^®(l + l|V4'||pa9i([;x(0,t)) + ll'^'t|lL“9i(;7x(0,d))^^ • ■‘^^) 


■ exp 


Concerning the power oft we note that S3,S2 > 1 and 


*2 *^1*3 \ 2g2W*3 _ 


2 — 


= S2P0S3 > pQ. Then 


(iT/mi follows 

(ii) Consider t > 3. Applying (|7.23l) with To = t — l,r=l,0 = 1/2 yields 

l|Vp(t)||Loo(t/) < c(l + / [ \Vp\''^dxdt] +C sup ||Vp||ioo(r). (7.26) 


Thanks to ^6.291) with s = S2, and (15.391) . and also noting that sup[i_3/4 K, 2 { ) < ^ 2 {t), we obtain 


INpimuo^iu) < C^[(l + \\Po\\LHU)f^^}^iitf^^^2{tf^^ (l + Gi(r)d- 


• exp 


{c(l + ||po||/)L(t;))/Cl(i)^^(l + 


IVTI 


L°^<n{Ux{t-2,t)) 


+ 11^ 


t||L“9i(f/x(t-2,t 




1*3/2 


+ CiC 2 {t) exp |c'^(l + ||po||l,“(!7))^®^l(i)(l + l|V'I'llL“9i([/x(t-2,t)) + ll'I't||l,“9i(;7xd-2,t)))^^|- 
Then p7.22p follows. 

Combining (I7.26P with Theorem 16.61 we have the following asymptotic estimates. 


□ 
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Theorem 7.3. (i) If A{a.) < oo then 
limsup||Vp(t)||ioc(c;) < 

t^OO 

■ exp [c'Ai{\ + limsup(||V'I'||j;^aiji([/xt_i,t) + (7.27) 

^ t^OO ^ 

(ii) If /3{a) < oo then there is T > b sueh that for all t > T, 

, 53/2 


W^Pimi^iU) < ClCiitr^^^lC 2 {tr^^(l+ [ Gi(r)dr) 

^ Jt- 3/2 ^ 


exp 


(c/x(i-3,*)))"=}• (7.28) 


Proof. Taking the limit superior as t 00 of (|7.26p . and using estimates (j6.38l) and (I5.40p give 


limsup||Vp(t)||ioc([;) < 


t^OO 


• exp 


ic'AYY + limsup(||V4'(t)||2,a9i([;x(t-i,t)) + ll'^'t(^)llL“'Ji(f/x(t-i,t))))^^| 

'' i—>-oo 


1 ^ 3/2 


+ CA 2 exp{C'.4^®(l + limsup(||V4;'(t)||iagi([/x(i-i,t)) + ll^i(^)llL“9i(;7x(t-i,t))))^^}- 

t—>-CO 


Then p7.27p follows. 

(ii) Combining p7.26p with estimate p6.4ip for s = S 2 , and with estimate (I5.41h we have 


l|Vp(t)|| 


L^(U) < 


c[ sup sup /C2(-)^^^fl+ [ Gi(r)dr) 

[t-3/4,t] ft-3/4,tl ^ Jt-7/4 ^ 


[t-3/4,t] 


t-7/4 


• exp{C' sup /Ci(-)^®(1 + II V4l'(t) 11^091 ([/x(t-ll/4,t)) + II 'I'i(0llL“'Jl(;7x (4-11/4,4)))^^} I 

[t-3/4,t] 


S 3/2 


+ c*(i + sup (ll'i'tll^oo + II^^IIl°° + ll^^^lli°°)) 

[t-7/4,t] 


•exp{C(l+,d(a)“ 2“+||T(a, Oil (1 + II^'^'(*)IIl“« (!7x (t-7/4,t)) + ll'^'t(^)llL“M(c/x(t-7/4,t)))'^^} 


for all t > T with some T > 3. This leads to P7.28I) . 


□ 


8 L'^-estimates for the time derivative 

We now estimate the L°°-norm of the pressure’s time derivative. Let p{x, t) be a solution of IBVP 
and let q = Pt- Then 

^ = V-{K{\Vp\)Vp)^. (8.1) 

Proposition 8.1. //Tq > 0, T > 0 and 9 G (0,1), then 

sup ||pOlL°°(t/) < CIA 2 (1 + (0r) ^) 2 ||pf ||L2((7x(ro+6»r/2,To+r)) + ^ 

[To+6»T,To+T] rxiro,To+TJ 


where si and A = X{To,T,6) are defined in Theorem 7.1. and constant C > 0 is independent ofT^, 
T, and 6. 
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Proof. Without loss of generality, assume Tq = 0. For k > maxpx[o,T] \Pt\ = ™axpx[o,r] IV’tl i let 
gi{fc) = max{g — A;,0} and Sk{t) = {x & U : q{x,t) > k}, and Xk{x,t) be the characteristic function 
of set {{x,t) G t/ X (0, r) : q{x,t) > k}. On Sk{t), we have {Vp)t = Vg = Vq^^\ 

Let C = C{t) be the cut-off function satisfying ^(0) = 0. We will use test function noting 

that Multiplying (I8.1|) by and integrating the resultant on U, we get 

lliJu = X " JjK{\Vp\))tVp • Vq^^\^dx - K{\Vp\){Vp)t • Vq^^\‘^dx. 

Note that (Vp)t • = |V(g'l^lC)P- Taking into account (12.121) . 

\iKi\Vp\))tVp • Vq^^k^l = \K'i\Vp\)\^-^^^^\Vp ■ < aK{\Vp\)\Vq\\Vq^^^C\ 

Note that |V(?|It follows that 

^ ^ \q^^\\^dx + (1 - «) X Ki\Vp\)\V{q^^\)\^dx < 2 ^ \q^^^C\Ct\dx. 

Integrating this inequality from 0 to T, we obtain 

max [ \q^’^^C\‘^dx+ [ K{\Vp\)\V{q^^\)\‘^dxdt < C [ \q^’^^\‘^C\Ct\dxdt. 

[0,T] Ju Jo Ju Jo Ju 

The last inequality uses the fact that function K{-) is bounded above. Applying Lemma [2.21 to 
with W = K{\Vp\) and Q = Vo defined by (I7.12p . we have 

lk^"^CIlL^o(Q^)<CAV-o./max / \qik)c\Hx+ r j K{\Vp\)\V{q^>^'>C)\^dxdt] ^ . 

I [OT] Ju Jo Ju J 

and hence, 

llq^^kh-'o(Q,) £lq(^^l%lCdxdt^^^^ (8.3) 

This is similar to inequality (I7.14p . Then by following arguments of Theorem 17.11 applied for 
instead of Um\ we obtain the estimate (18.2p . We omit the details. □ 

In the following estimates of pt we use the same notation as in sections [6] and 0 particularly, 
(I6.24p - (l6.27p . P6.39I) . ^6.401) . p7.17p - (|7.20p . and also define new numbers 

K4 = 1 -h Si -h S2 Ki(s3 - 1), K 5 = (S 2 Z 3 + 1)(S3 - 1) + 1, 

K6 = (S2Z3 + a/2)(s3 - 1) -h a/2 = ^3(53 - 1) -h a/ 2 . 

Theorem 8.2. (i) If 0 < t < 3 then 


\\Ptit)\\L--{u) <Ct + || po || l “) ^ (1 +j^H{\Vpo{x)\)dx^ ' 

■ G3{T)dT)y"^^ 


•exp 


[0,i] 


(8.4) 
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(ii) If t > 3 then 


wptimL^iu) < c(i + wpoh^r • 

{ C {l + \\^'^\\L°^n (f/x(t-3,t)) + ll'I'4|lL“M(f/x(t-3,t)))^^|+ SUp 


Proof, (i) Let t G (0, 3]. Applying (|8.2|) with Tq = 0, T = t and 0 = 1/2 gives 

si l + sj 

\\Pt{t)\\L-° < CX^t ^\\pt\\L^Ux{0,t)) +SUp||V’t||L“, 

[0,t] 

where A is defined by (j7.8p . By Young’s inequality, we have 


(8.5) 


( 8 . 6 ) 


^si/2 ^ 


t/4 JU 


{1+ \\Ip\)‘^-‘‘o dxdT] < 


t/4 Ju 


(1 + {Vpiy^dxdT 


where z^i = (2 — so)si/so = S 3 — 1- Combining with (j6.28p for s = S 2 yields 

+ J^Giiydr'^ 

■ exp ^C't + \\po\\la(^^'^y^)Cl{ty^{l + ||V^||x,a9l(;7x(0,t)) + ll^t||L“9l(f/x(0,t)))^^| 


Combining this with (18.61) and (13.91) gives 


pt 


exp{c't-G^yi + \\po\\L‘^iu)y^)cytyyi + ||v^ 

[JjH{\Vpo{x)\) +pl{x)]dx + GsiydrP 


0 

Il,“9i(t7x(0,t)) + ll^t||L“«(f/x(0,t; 
t r s 1/2 


1 / 1/2 



0 Ju 


»)"} 


9 \ 

^tix,T)\ dxdr) + sup 
^ [o,t] 


Therefore, 


||pt(t)||L- + 


(^ 22 : 3 + 1 ) 1 / 1+1 


iCi{ty^^^‘'^iC2{ty^''yi + 

exp 


^ G3(r)(ir)) 


1 + 

(i/l+l )/2 


'U 


Hi\Vpoix)\)dx 


1/2 


ict ^'^^^(1 + ||po||l“(;7))^^^i(^)^^( 1 + l|V'I'l|L“'?i(c/x(o,t)) + ll^t||L“9i([/x(o,t)))^^| + sup ||V’t||i,°°- 

J [0,t] 

Thus, we obtain (18.4h . 

(ii) Now consider t > 3. Applying (18.21) with Tq = t — 1, T = 1 and 0 = 1/2, we obtain 


||Pt(t)||L- < CX 2 ||pt||l,2([;x(t-3/4,t)) + sup llV'tlU-, 


(8.7) 
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where ^ ft- 3 / 4 fu(^ dxdtj . Then by Young’s inequality and (I6.29P for s = S 2 , 

^ 1/2 


1 / 1/2 


J (^i + \s/p\y^dxdT 

<c7(i + iipoiiL-r"^/ci(tr"i/c2(tr"^ • (i+J' ^Gi{T)d 

■ exp |c''(l + ||po||l,“(C/))^®^l(i)^®(l + l|V'I'||£,a9i([7x(4_3,t)) + ||'I't||l,“9i((7x(t-3,t)))^^|- 
Combining this with (18.7p and (I3.1ip gives 

ui/2 




\\pt{t)h^(U) < C(1 + • (1 + l_^GyT)dT 

■ exp ^C'{1 + ||po||L“(i7))^®^l(i)^®(l + l|V'I'||^a9i([7x(t_3,t)) + ||'I't||L“9i (17x0-3,0))^^} 

•( 1 + f \po{x)\°'dx + [EnvA{a,t)]°^ + [ G 3 (r)(ir+ f f l^tl’^dxdr) ^ + sup 
^ Jv Jt -2 Jt -2 Ju ' [i-l,i] 

Using the facts Gi{t) < G^^t) and Jjj \'i>t{x, t)\‘^dx < G 3 {t), we infer that 
\\ptmL^iu) <c{i + \\Po\\L‘^r'^^^‘^G)cytr''i+-G)cytr^^ 

and obtain (j8.5p . The proof is complete. 

For large time or asymptotic estimates, we have the following. 

Theorem 8.3. (i) If A{a) < 00 then 

limsup||pt(t)||ioo(^) < 

^ vco 

) I G'Ay (1 + lim sup (11VT11 ^091 (f; X (t-1 ,i)) + 11 (t) 11 L“91 ([/ X (7-1 ,t))) )^^ I + lim sup 11 112.“ (r), 

7 —>00 7 i —>-00 


□ 


t^OO 

• exp 


( 8 . 8 ) 


where ^2 = 1 + limsupj_^(x) //_i G'i{T)dT. 

(ii) If /3{a) < 00 then there is T > 0 sueh that for all t > T, 


\ss/2 


/ /*^ \ * 
lb(i)llL-(C/) / GyT)dT\ 

(8.9) 

•exp|c"Yi(t)^’^(l + ||V4'||i/«9i(;7xO-3,t)) + \\^t{t)\\L‘-n{ux{t-3,t))y^] 


Proof. Again, let z^i = S 3 — 1. 
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(i) Taking the limit superior of ()8.7I) as t —>■ oo, and using (|6.38p with s = S 2 , and (13.131) . we 
obtain 


limsup ||pt(t)||Loo([/) <C 


t^OO 


• exp 


ic'A^iil + hmsup(||VT||2,a9i(t/x(t-i,t)) + \\^tit)\\L°‘ii(ux{t-i,t)))y^\ 

t—>0O 


11/1/2 


Af + lim sup 

i-l-oo Jt-l 


/ G 3 (r)dr + limsup / / 

J t —1 t^oo J t —1 J U 


1/2 


Tt(x,r)pdxdTj + limsup \\yt\\L°°{r)- 


t—>-oo 


Simple manipulations yield 


limsup ||pt(t)||ioo(( 7 ) < CA 


t^OO 

• exp 


^ ^ ^ Hm g^p / G^{T)dT 

t^oo Jt-l 


1/2 


I C'Ay (1 + lim sup (11VT117,091 (t; X (t -1 ,i)) + 11 (i) 11 ([/ X ( 1-1 ,t)))) ^^ I + lim sup 1111(r). 

t—>-oo ^ t—>-oo 

( 8 . 10 ) 


t—>-oo 


Since Qi < G 2 , we have + limsup^^oo//_i ^ Thus, 

inequality (18.8h follows (I8.10h . 

(ii) Combining inequality (18.7p with estimate (I6.41|) for s = S 2 and estimate (j3.15p . we obtain 
for large t that 




exp 


|c''^l(t)^3(l _|_ lIVTllT^oqij-T/x/t-S,!)) + ll'l'i(^)llL“«(;7x(l-3,0))^^} 


1/1/2 


1 +/3(q;)“- 2“ + sup /4(a, r — 1) + / G 3 (r)dr + 




n -2 


lt-2 Ju 


|Ti(x, T)p(ixdr 


1/2 


+ max hAt L°°(n- 

Using similar calculations as in part (i) and the fact Gi < G 3 , we obtain (|8.9|) . 


□ 
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